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Remote Minehunting 





An image of the target field using the high frequency synthetic aperture sonar. The sonar is 25 feet 
above the bottom and the water depth is 98 feet. The maximum range displayed in the image is 
approximately 250 feet. In the near ranges, a number of cylindrical targets of varying sizes and 


orientations are clearly visible, as are their shadows. Also evident are the images of two spheres which 
appear as a pair of bright dots. 
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Foreword 


Douglas Todoroff, Guest Editor 
Office of Naval Research 


In the post-Cold War environment, naval warfare pri- 
orities have shifted from open/blue-water-basin operations 
to regional, shallow water Third World conflicts as the ma- 
jor technology challenges facing naval operations in the 21st 
century. This increased importance of responding to littoral 
conflicts has highlighted the demanding requirements for 
developing effective sea-mine countermeasures in shallow 
and very shallow water. Thus, it is appropriate at this time 
that the Naval Research Reviews should devote a theme is- 
sue to the critical problem of countering the mine threat. In 
this issue, we focus on recent technology developments in 
the area of mine detection, classification and identification. 

The recent interest in mine countermeasures is a result 
of the proliferation of mines and the threat they pose to na- 
val operations. Mines can be easily obtained and used by 
any nation. Sophisticated mines can now be purchased com- 
mercially, and crude copies of sophisticated mines can be 
manufactured in-house, even by those nations with limited 
industrial capabilities. As a result, mines have become a 
weapon of choice for Third World countries. Mines have 
been used repeatedly in the last quarter century to threaten 
commercial shipping, to cover aggressive actions, and to at- 
tack naval vessels. In recent years, USS Roberts, USS Tri- 
poli, and USS Princeton have all been victims of naval min- 
ing. 

The diversity of mines available to regional powers 
makes shallow water mine countermeasures perhaps the most 
technology challenging of all naval mission requirements. 
In deep water, tethered mines constitute the main threat to 
naval operations. In shallow water , however, the mine threat 
consists of bottom, volume, and floating mines. Here shal- 
low water refers to depths less than 100 meters. A very shal- 
low water regime has been defined as the coastal region with 
water depths from 11 to 3 meters. In this regime, the threat 
consists primarily of bottom-influence sea mines, smaller 
anti-invasion mines, moored mines, and floating mines. 
Additionally, the likelihood of mine burial increases as the 
water depth decreases. 

As a result, there is no single, simple solution to the 
mine detection problem. In particular, it is generally recog- 
nized that pure acoustic approaches to the detection of un- 
derwater sea mines cannot be as effective as multidisiplinary 
techniques. In this issue, we include recent mine detection 
technology advances in both acoustics and non-acoustics. 

The first article is by researchers at Coastal Systems 
Station (NSWC) and Northrup Grumman regarding recent 
developments in synthetic aperture sonar for mine detection 
and classification. Synthetic aperture sonar, like synthetic 
aperture radar, involves the processing of data from a small 
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physical aperture to provide high resolution imagery. In this 
article, a description of a novel motion compensation tech- 
nique is presented which is critical to the application of syn- 
thetic aperture sonar technology for minehunting. 

In the second article, Dobeck, et.al. describe an ap- 
proach to automated detection and classification of sea mines 
in sonar imagery which has performance similar to that of 
an expert sonar operator. Their approach utilizes an improved 
detection density algorithm, a feature extraction technique 
that uses a stepwise feature selection strategy, and two inde- 
pendent classifiers. In this article, the authors demonstrate 
significant improvement in overall classification by 
“ANDING?” classifiers which use complementary statistical, 
mathematical, and geometrical constructs to describe class 
boundaries in feature space. 

The third article address the identification of proud, 
partially buried, and moored mines in shallow and very shal- 
low water. Mine identification is a critical phase of 
minehunting in which mines are discriminated from mine 
like contacts prior to neutralization. In this article, Mike 
Strand describes a comparison of laser range-gated and la- 
ser line scan technologies for mine identification. 

The final article describes recent advances in the mag- 
netic detection and classification of sea mines and of 
unexploded ordnance. In this article, Ted Clem describes 
the development of a 5-channel tensor magnetic gradiom- 
eter that provides significant improvements in classifica- 
tion and localization over that offered by the commonly used 
single-channel total field magnetometer. While the current 
configuration of this sensor features bulk niobium low criti- 
cal temperature superconducting components, recent ad- 
vances in material research will allow development of a high 
temperature superconducting gradiometer concept using liq- 
uid nitrogen refrigerant to reduce package size and cryogenic 
support requirements. 

The technologies described in this issue of Naval Re- 
search Reviews represent significant advances in the detec- 
tion, classification, and identification of sea mines. The syn- 
thetic aperture sonar, automated processing, and electro-op- 
tic identification technologies will be demonstrated to the 
operational user as part of the Joint Countermine Advanced 
Concept Technology Demonstration (ACTD) in fiscal year 
1998. 





High Frequency/Low 
Frequency Synthetic 
Aperture Sonar 


Gary S. Sammelmann, Jose E. Fernandez, James T. Christoff, Leonid Vaizer, & John D. Lathrop 
Coastal Systems Station, Dahlgren Division, Naval Surface Warfare Center, Panama City FL 


Robert W. Sheriff & Thomas C. Montgomery 


Northrop Grumman, Oceanic Division, Annapolis MD 


Abstract 


The HF/LF SAS (high frequency/low frequency SAS) 
is a high resolution SAS developed by COASTSYSTA and 
Northrop Grumman for the shallow water (SW) and very 
shallow water (VSW) regimes. This sonar suite has recently 
been delivered to COASTSYSTA and it is currently under- 
going field testing. This article describes this sonar and the 
type of resolution and acoustical images which are expected 
from this sonar. The application of this sonar to the SW and 
VSW regimes required the development of a novel method 
of motion compensation. A description of this method and 
the type of accuracy’s expected from this technique are pre- 
sented. Finally, a look at future broad band systems and their 
predicted performance is presented. 


1. 


Introduction 


The HF/LF SAS ( high frequency / low frequency syn- 
thetic aperture sonar ) is a high resolution synthetic aperture 
sonar developed by COASTSYSTA and Northrop Grumman 
to detect, classify, and identify mine-like targets in the shal- 
low water (SW) and very shallow water (VSW) regimes. 
The LFSAS is expected to be able to detect and classify vol- 
ume, proud, and buried mines with its 3x3” resolution. The 
HFSAS with its 1”x1” resolution is expected to have an iden- 
tification capability approaching that of some optical sys- 
tems. 

The HF/LF SAS was delivered to COASTSYSTA in 
July of 1996, and it is currently undergoing field testing. This 
article describes the design and the type of performance ex- 
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pected from the HF/LF SAS. This article also describes the 
novel motion compensation technique developed by 
COASTSYSTA and Northrop Grumman for this sensor. 

Section 2 of this article describes the design of the HF/ 
LF SAS, and Section 3 describes the motion compensation 
algorithms. Section 4 describes the current performance of 
the HF/LF SAS 


2. 


Hardware Design 


The main design objective of HF/LF SAS was to de- 
sign a small, light weight, high resolution acoustic sensor 
for operation in the SW and VSW regimes which was com- 
patible with operation in a UUV (unmanned underwater ve- 
hicle). Compatibility requirements with MK 48 based UUVs 
constrained the dimensions of this sensor to a cylindrical 
body 21 inches in diameter in a section no longer than 36 
inches. The requirement that the sensor be able to detect, 
classify, and identify mine-like targets and fit within the above 
section can only be satisfied by a SAS. The requirement that 
it be able to detect and classify buried targets, and identify 
volume and proud targets suggested a dual frequency ap- 
proach consisting of a low frequency SAS similar to 
MADOM to detect and classify buried targets, and a high 
frequency SAS for identification of volume and proud tar- 
gets. The horizontal and vertical beam widths of the sensor 
and the depression angle of the sensor have been optimized 


to produce the greatest signal-to-noise ratio (SNR) and cov- 
erage in the shallow and very shallow water environment 
using PC SWAT. The length of the receiver was chosen to 
attain the desired maximum range at a given platform veloc- 
ity. 

Past experience with MADOM showed that a LFSAS 
with a central frequency of 20 kHz and band width of 10 
kHz was capable of detecting and classifying buried targets. 
As a proof of concept, the original MADOM configuration 
was modified to have a central frequency of 20 kHz with a 
10 kHz band width and 6 inch elements. This configuration 
gave MADOM the capability of achieving 3”x3” resolution. 
In December 1992, a trial sea test of this configuration dem- 
onstrated the concept. The concept was further demonstrated 
in September 1995 at the MUDSS (Mobile Underwater Sur- 
vey System) Feasibility Demonstration, where this modified 
version of MADOM successfully detected and classified tar- 
gets ranging in size from a 60 millimeter howitzer shell to a 
2000 Ib. bomb. During this demonstration, the HPSS ( High 
Performance Sidescan Sonar ) was employed in a SAS con- 
figuration as a HFSAS. However, due to problems with the 
filters, this sonar was only able to achieve 5 cm x 15 cm 
resolution. 

The major obstacle inhibiting performance of the LF/ 
HF SAS was the motion of the sensor away from a straight 
line trajectory. The method of overcoming this problem was 
to develop an adaptive motion compensation technique which 
determines the ping to ping motion of the vehicle based on 





Figure 1 
Photograph of the completed system. 
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the correlation between elements of neighboring pings. The 
method chosen for this motion compensation was a tech- 
nique whereby the last 2 phase centers of the array overlap 
with the first 2 phase centers of the next ping. The time de- 
lay estimates between these overlapping pairs of phase cen- 
ters allows one to determine the translation, pitch, and yaw 
motion of the array between adjacent pings. For this pur- 
pose, the array was designed with an extra pair of elements 
in order to facilitate the necessary 2 phase center overlap 
between adjacent pings and reduce the magnitude of the grat- 
ing lobes of the array. 

To facilitate the motion compensation estimates of the 
array, a Motion Measurement Package (MMP) consisting of 
a Doppler sonar, and an Inertial Navigation System (INS) 
was included. The INS allows one to determine an approxi- 
mate value of the motion of the sensor. However, the accu- 
racy of INS is insufficient to achieve the desired resolution 
of the sensor. Hence, the need for an adaptive motion com- 
pensation technique, which has the potential of estimating 
the motion of the vehicle to within a thousandth of a wave- 
length, is necessary. However, the INS provides the operator 
with a baseline to compare the estimated motion with the 
motion measured by the INS. The Doppler sonar provides 
accurate velocity measurements. By slaving the ping repeti- 
tion rate of the projector to the Doppler unit, one can reduce 
the effects of time varying velocity of the array along its 
trajectory, by requiring that the projector fire after the array 
has moved a predefined distance along track. 


The projector of the LFSAS is a tonpilz design. The 
vertical height of the projector can be configured as either a 
one or two wavelength high projector by changing the jumper 
leads on the matching network assembly. The projector can 
produce a high fidelity 2 cycle CW pulse with a 10 kHz band 
width centered at 20 kHz. 

The projector of the HFSAS is a standard half wave- 
length ceramic bar design. The projector consists of 10 ele- 
ments in a linear array, which allows one to vary the center 
of the projector at will. The projector can produce a 6 cycle 
CW pulse with a 30 kHz bandwidth centered at 180 kHz. 

The receive elements of the LFSAS and HFSAS have 
been designed using a PiezoRubber (PZR) as the active ma- 
terial. The LFSAS contains 14 elements with an element 
spacing of 1.5 inches. Only 12 of these 14 elements are used 
to form the synthetic beam. The remaining 2 are used to 
achieve the desired phase center overlap between adjacent 
pings. The HFSAS contains 11 elements with an element 
spacing of 2 inches. Only 9 of these 11 elements is used to 
form the synthetic beam. The remaining 2 are used to obtain 
the desired phase center overlap. 

The number of elements in the LFSAS and HFSAS is 
determined by the maximum range (40 meters) of the sen- 
sor, and the maximum platform velocity (8 knots). The dis- 
tance traveled between pings is 0.2286 meters or 9 inches. 
The LFSAS forms 3 beams with a spacing of 3 inches per 
ping, and the HFSAS forms 9 beams with a spacing of 1 
inch per ping. 





Figure 2 
Exploded view of the 21 inch section. 
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A study of the predicted signal-to-noise ratio of the 
LFSAS and HFSAS as a function of depression angle showed 
that the signal to noise ratio for a proud target was maxi- 
mized and the blind spot under the vehicle was minimized 
for a depression angle of 10 degrees. 

Figure | shows a photograph of the completed sys- 
tem. Figure 2 shows an exploded view of the 21 inch sec- 
uon. 


3. 


Motion Compensation 


The method of motion compensation developed by 
Northrop Grumman and COASTSYSTA is an adaptive ap- 
proach which utilizes the cross-correlation between elements 
of adjacent pings to estimate the time delay. This approach 
requires that the phase centers of the last 2 elements of a 
ping and the first 2 elements of the next ping overlap. Two 
overlapping elements are required to determine the ping-to- 
ping time delay and its rate of change. 


The concept of using overlapped displaced phase cen- 
ters (DPC) was first disclosed by R.S. Raven [1]. Subse- 
quent work at Northrop Grumman and COASTSYSTA [2, 3] 
has demonstrated the ability of this technique to compensate 
for ping-to-ping translations and rotations. 

The central #dea to this motion compensation technique 
is the assumption that the signal from a pair of overlapping 
phase centers should be nearly identical, except for small 
deviations due to the motion of the platform. This assump- 
tion is based upon the notion that both phase centers are illu- 
minating approximately the same footprint on the bottom, 
and have approximately the same travel time to any point on 
the bottom. However, contamination of the signal due to 
multipathing, and surface reverberation in particular, at the 
larger ranges tends to decrease the correlation between the 
two signals, since multipathing adds a time varying signal to 
the signal at the phase centers. Thus, one must compensate 
for the effect of multipathing at the larger ranges. 

The second key assumption to this motion compensa- 
tion technique is that the signals from the overlapping phase 





Figure 3 

Point response function of the HFSAS with no motion in the 
three cases: no motion compensation; one element overlap; 
and two element overlap. 


Figure 4 

Point response function of the HFSAS with with a constant 
10 degree crab in the three cases: no motion compensation; 
one element overlap; and two element overlap. 
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centers have a high correlation, and one can estimate the 
time delay between the two signals due to the motion of the 
platform by estimating the time delay in the cross correla- 
tion between these two signals. By using two overlapping 
phase centers, one can estimate the time delay and the rate 
of change of the time delay between two adjacent pings due 
to the motion of the platform. 

Figure 3 shows the point response function of the 
HFSAS with no motion in the three cases: no motion com- 
pensation; one element overlap; and two element overlap. In 
all three cases, the point response function is close to its theo- 
retical value. 

Figure 4 shows the point response function with a con- 
stant 10 degree crab in the three cases: no motion compen- 
sation; one element overlap; and two element overlap. The 
case with no motion compensation is greatly distorted by 
the motion, whereas the two cases with motion compensa- 
tion are similar to those without motion, indicating the mo- 
tion compensation technique has properly removed the dis- 
tortion of the synthetic beam pattern caused by the constant 
crab. This case shows that the case of one element overlap is 
sufficient to remove static time delays due to motion. 





Figure 5 
Point response function in the case of a sinusoidal yaw. 
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Figure 5 shows the point response function in the case 
of a sinusoidal yaw. The case with no motion compensation 
is greatly distorted by the yaw. The case with a one element 
overlap is somewhat distorted due to the time varying nature 
of the yaw. 

The case with a two element overlap is similar to the 
case with no motion indicating that the motion compensa- 
tion technique with a two element overlap has correctly re- 
moved the distortion to the image due to the time varying 
yaw. 


4. Performance 


Figure 6 depicts an image of the target field during the 
sea tests of the HFSAS in the shallow water field. 


The picture on the inside front cover of this issue de- 
picts an image of the target field using the HFSAS. The 
image of 2 spheres and the resolution panel are clearly in 
evidence in this figure. The image of each sphere consists 
of 2 returns. The specular return, and a second return which 
may be an elastic return due to penetration of the fluid filled 
spheres. The legs of the resolution panel are clearly in evi- 
dence in this image, as are the doughnut shaped protrusions 
lying on the resolution panel. 





Figure 6 
Image of the target field using the HFSAS. 
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5. Conclusions 


The recent completion of the HFSAS and LFSAS and 
the September 1996 sea test of this system offers new capa- 
bilities in mine detection, classification, and identification 
in shallow and very shallow water. The motion compensa- 
tion technique described in this article holds the potential of 
predicting the time varying motion of the HFSAS and LFSAS 
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to a high degree of precision ( a thousandth of a wavelength). 
During the next few months further refinements of this mo- 
tion compensation technique will be worked out to refine 
the process and improve the quality of the HFSAS and 
LFSAS images. 
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Abstract 


An advanced capability for automated detection and 
classification of sea mines in sonar imagery has been devel- 
oped. The Advanced Mine Detection and Classification 
(AMDAC) algorithm consists of an improved detection den- 
sity algorithm, a classification feature extractor that uses a 
stepwise feature selection strategy, a k-nearest neighbor 
attractor-based neural network (KNN) classifier, and an op- 
timal discriminatory filter classifier (ODFC). The detection 
stage uses a nonlinear matched filter to identify mine-size 
regions in the sonar image that closely match a mine’s sig- 
nature. For each detected mine-like region, the feature ex- 
tractor calculates a large set of candidate classification fea- 
tures. A stepwise feature selection process then determines 
the subset features that optimizes probability of detection 
and probability of classification (PdPc) for each of the clas- 
sifiers while minimizing false alarms. 

The AMDAC has been tested using 335 sonar images 
from three different sonar systems: a synthetic aperture so- 
nar (SAS with 255 images), a side-scan sonar (SSS1 with 60 
images) and another side-scan sonar (SSS2 with 60 images). 
The AMDAC’s performance is: 90% PdPc with 0.42 false 
alarms per image for the SSS1; 92% PdPc with 0.64 false 
alarms per image for the SAS; and 91% PdPc with 0.12 
false alarms per image for the SSS2. For these data bases 
the algorithm’s performance is as good or slightly better than 
that of an expert sonar operator. 

This work was funded by the Office of Naval Research, 
ONR 32ITS, through the 6.2 Mine Countermeasures pro- 
gram element. 


S; 


Introduction 


This paper builds on our earlier work [1] and describes 
a new algorithm for the automated detection and classifica- 
tion (D/C) of sea mines in sonar imagery. Three sonar im- 
age data bases are used to evaluate the new D/C algorithm. 
The algorithm is referred to as the Advanced Mine Detec- 
tion and Classification (AMDAC) algorithm. It is comprised 
of four stages: (1) Image Enhancement, (2) Detection, (3) 
Feature Extraction and Optimal Feature Selection, and (4) 
Classification (see Figure 1). 

For the last seven years CSS has used three side-look- 
ing sonar data bases to evaluate mine detection and classifi- 
cation (D/C) algorithms developed in house, by industry and 
by academia [2], [3], [4], [5]. The data bases consist of a 
255 image set from a synthetic aperture sonar (SAS), a 60 
image set from a side-scan sonar (SSS1), and a 60 image set 
from another side-scan sonar (SSS2). 

These data bases were selected because they demon- 
strate a variety of problematic issues and signal processing 
challenges. Mine threats in these data bases are bottom 
mines; they provide a significantly greater challenge to de- 
tect and classify than volume mines. Results for probability 
of detection and classification (PdPc) and false alarm rates 
referred to in this report are for single pass (one look) sonar 
operation. Single pass operation is of interest because of the 
implied higher search rate. Performance of an expert sonar 
operator was available for two of these databases, thus pro- 
viding a benchmark for automated algorithm performance. 
Because these three data bases are somewhat small and do 
not represent all interesting scenarios, one must be cautious 
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Figure 1 
Overview of AMDAC Algorithm. 
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in generalizing any conclusions. However, the prohibitive 
cost of collecting and processing large data bases makes 
utilization of small data bases the typical practice when evalu- 
ating and comparing automated D/C algorithms. Real so- 
nar image data is preferred over simulated sonar data be- 
cause sonar simulations are expensive and do not capture all 
the critical dynamics associated with actual sonar images. 

The following sections, 2.1 through 2.6, describe the 
four stages of the AMDAC. 


2. AMDAC Algorithm 


2.1 Image Enhancement 


The purpose of the image enhancement stage is to pre- 
condition the image so that the subsequent detection and clas- 
sification stages are robust to variations in background level. 
This is accomplished by normalizing the background 
throughout the image to a constant level so that highlight 
and shadow levels are consistent and clearly stand out. 

The principal reasons for variations of bottom bright- 
ness are (1) inadequate range-varying gains that are applied 
to the sonar return to compensate for the range-dependent 
reduction in signal strength cause by spherical spreading of 
acoustic energy, (2) the strong sonar returns from highly re- 
flective bottom regions (e.g., gravel), and (3) the weak sonar 
returns from highly absorbent regions (e.g., mud). Thus, pre- 
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conditioning reduces some of the background invariance and 
reinstates some degree of robustness to bottom variation. 
For the three sonar data bases studied in this article 
the inadequate range-varying gain problem, as described 
above, was the dominate factor; so a simple range normal- 
ization was found to be adequate. This was done as follows. 


I (i,j) = raw 8-bit sonar image (0 < I (i,j) < 255 ) 
i = range index (Ilsisn) 
j = cross-range index (1<Sj<m) 
I (i,j) = normalized image (0 <1 (i,j) <4) 
= MIN (4, (i,j)/(i) ) 


where 


b(i) = mean level of I at range i 


Note the following facts about I : 


(1) The mean background pixel intensity is essentially unity. 
(2) The dynamic range for highlight pixel intensity is from | 
to 4. 





(3) The dynamic range for shadow pixel intensity is from | 
to 0. 


Experience of expert side-looking sonar operators has 
shown that the above scaling procedure is a good way to 
distribute the dynamic range between highlight and shadow. 

In subsequent discussions we will make use of these 
three facts in setting parameter values. As an aside, it has 
been found that images from other environments may con- 
tain many large irregular regions, some of which are highly 
reflective while others are highly absorbent. In these cases 
the simple range normalization does not achieve the goal of 
making the background level uniformly near unity. Other 
computationally more intensive methods are needed that 
normalize the image based on local statistics within the im- 
age. As an example, low-pass median filtering has performed 
quite well for us. 


2.2 Detection Stage 


The purpose of the detection stage is to scan the entire 
image and identify candidate mine-like regions that will be 
more thoroughly analyzed by the subsequent classification 
stages. This is the most computationally intensive stage be- 
Cause a mine-size region surrounding each image pixel must 
be evaluated. Therefore, the goal is to keep the computa- 
tions involved with each mine-size region small. The pur- 
pose of this computation is primarily not to determine mine- 
likeness but rather to screen out the non mine-like regions in 
the vast majority of the image and thus reduce the amount of 
data that must be processed by the computationally more 
intense classification stages. If the detection stage can re- 
duce a typical image to about 10 detections (10 being very 
generous), the computational requirements of the classifica- 
tion stages are insignificant when compared to the detection 
screening. 

The detection stage of the AMDAC is an improvement 
on the one reported in our earlier paper [1]. The detection 
stage operates on the normalized image. The algorithm di- 
vides the image into regions along range to account for the 
variability of the background and mine signature as a func- 
tion of range. For this study the image was divided into 
three regions. Next, the image is convolved with a nonlin- 
ear matched filter, a different matched filter for each of the 
range regions. The procedure sets negative values from the 
match filter to zero. For each range region the resulting 
matched-filtered image is normalized by removing the 
region’s mean and dividing by the region’s standard devia- 
tion. Next, the procedure scans a target-size window over 
the normalized matched-filter image. A detection is declared 
by counting the total number of pixels within the window 
that exceed a specified “amplitude-detection threshold”; a 
pixel count above the specified “detection-count threshold” 
is declared a target. 


2.3 Nonlinear Matched Filter 


The nonlinear matched filter is the work horse of the 
detection stage [6]. The matched-filter mask contains four 
distinct regions: pre-target, highlight, dead zone and shadow/ 
post-target ; see Figure 2. It is defined as, 


I, (i,j) = y ¥ e(h(k,1),1, (1 +k, j+1)) 


=-N, l=-M, 


h(k,1) (I-1) for h(k,l) corresponding to 
shadow, highlight, and dead zone regions 


h(k,]) lI-1! for h(k,]) corresponding to pre- 
target region and post-target regions 


Note that the nominal background level of unity is subtracted 
from I. 





Figure 2 
Target Signature Mask. 
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In each of the four regions, the matched filter coeffi- 
cients are constant and defined by, 


shadow region or post-target region: h(k,l) = 
highlight region: 

dead zone region: 

pre-target and post-target regions: 


1MS,(S, - 1) 
1(H,(H, - 1)) 
0 


= -1(TIT, - 1) 
where 


S,= area of shadow region in square pixels (see Tables 
1, 2 and 3) 


S, = reference shadow level = 0.75 for this study 


H, = area of highlight region in square pixels (see Tables 
1, 2 and 3) 
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H, = reference highlight level = 1.5 for this study 


T, = area of pre-target region in square pixels (see 
Tables 1, 2 and 3) 


T, = reference anomalous background level 


The following two anomalous levels used in this study 
result in the same h(k,]) value. 


2.0, anomalous high backgound level 
0.0, an anomalous low background level 


In this study the reference center h(0,0) corresponds 
to the center pixel of the highlight region. The following is a 
brief discussion of this matched filter. 

The nonlinearity is associated only with the pre-target 
and post-target regions. The intent of this part of the matched 
filter was to prevent the detection of objects with highlights 
or shadows that were greater than mine size. The filter coef- 
ficient for this region is negative; thus, the output of the fil- 
ter is decreased if the pre-target or post-target regions are 
dominated by low or high valued pixels. This region could 
have been selected to encircle the typical mine-size area. 
However, the large variations in the mine’s cross-range size 
and shadow length, together with considerations of the extra 


computational burden, led us to use only the area in front of 


the highlight region. Since the mines in our SAS image data 
base have no shadows, a post-target region was used in place 
of the shadow region. 

The selection of the coefficients, defined by the refer- 
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ence levels and the region’s area, was done to make the con- 
tribution from each region of the matched filter equal for 
typical shadow, highlight, and anomalous background lev- 
els. Specifically, the coefficients are constructed such that 
the contribution for each region is unity when evaluated at 
the reference level of that region. 

The optimal detection-amplitude and detection-count 
threshold pair is selected from the optimal receiver operator 
characteristics (ROC) table [6]. The ROC table contains 
optimal threshold pair that gives the least false alarm rate 
for a given Pd. A computationally fast method was devel- 
oped to generate the ROC table that requires only a single 
pass through the data base. 

The AMDAC uses a different filter mask for each mine 
type and each range region. For example, the SAS sonar has 
one mask for detecting mine type “A” in the first region, 
another mask for detecting mine type “A” in the second re- 
gion, and another mask for detecting mine type “A” in the 
third region. To detect a different mine type, three different 
masks would be used. Reference [1] discusses target mask 
design in detail. Tables | through 3 show the masks’ sizes 
for the three sonars and the three range regions. 


2.4 Feature Extraction 


In the feature extraction stage, a large set of classifica- 
tion features is computed from a target-size classification 
window centered at each detection location in both the nor- 
malized raw sonar image and the normalized matched-filtered 
image. The procedure calculates forty-five candidate features 
based on the size, shape and strength of the highlight and 
shadow, and histogram information of pixel intensity within 
the classification window. A complete list of features fol- 
lows: 

. W_pix: number of normalized image pixels in the 
window that exceed a threshold* 
number of matched-filter pixels in the win- 
dow that exceed a threshold 
maximum normalized image pixel intensity 
in the window 
maximum matched-filter pixel intensity in 
the window 

target strength computed from the normal- 

ized image image 
target strength computed from the matched- 
filter image 
length of major axis of an ellipse fit to high- 
light region of the normalized image 
width of minor axis of an ellipse fit to high- 
light region of the normalized image 
length of major axis of an ellipse fit to bright 
region of the matched-filter image 


2. Mfw_pix: 

. M_wpix: 

. Mmf_wpix: 
. T_str: 

. Mft_str: 

. Mx_eig: 

. Mn_eig: 


. Mfmx_eig: 





10. Mfmn_eig: width of minor axis of an ellipse fit to bright 
region of the matched-filter image 

11. S_len: shadow length 

12. S_str: shadow strength 

13. Mmf_pclu: maximum matched-filter intensity over the 
pixels in the detection cluster 

14. C_pix: number of pixels in the detection cluster 

15. C_im: number of detected clusters in the image 

16 - 25. Hist_org(i): umber of normalized pixels above 
threshold(i) in the window, i from 1 to 10 

26 - 35. Hist_mf(i): number of matched-filter pixels above 
threshold(i) in the window 

36 - 45. Hist_diff(i): number of matched-filter pixels above 
threshold(i) in the classification window 
minus the number of matched-filter pixels 
above threshold(i) in the region that locally 
surrounds the classification window 
* All thresholds are selected by the designer. 


2.5 Optimal Feature Selection 


When training classifiers with finite databases, Bellman 
discovered that the robustness of a classifier will collapse 
when the number of features becomes too large [7]. This is 
known as “the curse of dimensionality.” Therefore, an im- 
portant part of the overall classifier design is to select a small 
and robust set of classification features from the larger can- 
didate set. Evaluating all possible combinations of 45 fea- 
tures is computationally not feasible (there are 2* - 1 combi- 
nations, a number greater than 10'*). Therefore, a stepwise 
optimal selection process is used (select the best single fea- 
ture, next select best of the remaining 44 features to add to 
the first, next select to best of the remaining 43 features to 
add to the previous two, etc.) [8]. For a 45 feature set this 
involves evaluating 1035 feature combinations. 
Computationally efficient algorithms have been developed 
for this selection process. 

The two classifiers used in this algorithm were the k- 
nearest neighbor neural network (KNN) and the optimal dis- 
criminatory filter classifier (ODFC). Since two classifiers 
are used in this detection/classification algorithm, there are 
two selection algorithms; each one is tuned to the respective 
classifier. This is required because a set of features that is 
optimal for one classifier will, in general, be far from opti- 
mal for the other classifier. This usually results when each 
classifier uses a different mathematical/geometrical struc- 
ture to partition feature space into class regions. A set of 
features amenable to one partitioning scheme may not sepa- 
rate well for a different partitioning scheme. 

Briefly, two feature selection algorithms were devel- 
oped for the KNN and ODFC, respectively. By optimal is 
met the following. Bellman’s curse of dimensionality im- 
plies that, for a finite training set, the detection and classifi- 
cation performance (as judged by the validation/test set) does 


not continually improve as more features are added (as di- 
mensionality is increased) to the classification process. With 
a finite data set it is possible only to determine a subset of 
features that will give robust performance. Optimal selec- 
tion in this context means finding the subset of features from 
a much larger set that gives the best performance for both 
the finite training and validation (test) sets for a specific clas- 
sifier (in our case, either KNN or ODFC). Specifically, the 
best subset is the subset of features which when fed into the 
classifier minimizes a risk metric. The risk metric that is used 
is defined as the maximum (between the training and valida- 
tion data sets) of the weighted sum of the number of missed 
mines and number of false calls. As stated above, for a typi- 
cal candidate feature set of size 45, there would be 2% - | 
combinations of subsets to evaluate in order to determine a 
global optimum. Because this is too large of a number to 
evaluate, a Forward Stepwise Optimal Selection Process 
(FSOSP) is used instead of a Globally Optimum Selection 
process. The FSOSP selects as the first feature the one that 
gives the best classifier performance of any single feature. 
The second feature selected is the one that performs best 
when used with the first. As it name implies, the Forward 
Optimal Stepwise Selection Process continues at each stage 
by adding a new feature to the subset determined in the pre- 
vious stage. This stepwise process is obviously not globally 
optimal, but has worked extremely well in this optimization 
application. Our research suggests that it is important that 
the classifier itself be used in the selection process. Other 
selection techniques (e.g., principal component analysis or 
multivariate normal models) are often used because they are 
computationally tractable. But they are usually not based 
on the same mathematical, statistical, or geometrical struc- 
ture as the classifier and, therefore, do not select the best 
features for that classifier. Keeping in mind Bellman’s curse 
of dimensionality, our research suggests that classifiers that 
train fast, for which the best subset of features can be opti- 
mally selected, will perform significantly better than more 
sophisticated classifiers for which the best subsets cannot be 
found because training is too computationally expensive. 

The selection process has been further improved by 
adding a Backward Stepwise Optimal Selection Process 
(BSOSP) analogous to model order reduction used in statis- 
tics [9]. Because the candidate features are not orthogonal, 
adding a new feature to the current optimal feature subset 
can in fact decrease the performance metric; this is consis- 
tent with Bellman’s curse of dimensionality. The augmented 
procedure is to first use the FSOSP to determine the initial 
feature set as previously described. Then, the BSOSP deter- 
mines which features to remove. Application of the FSOSP 
and BSOSP can be repeated until performance improvement 
stops. As an example, using the SSS2 data base, the num- 
ber of false alarms per image decreased from 0.28 (using 
only FSOSP) to 0.12 (using FSOSP followed by BSOSP) 
for the same PdPc of 91%. 


Three/1997 





2.6 Classification 


The KNN and ODFC classifiers were chosen because 
they process the data using significantly different mathemati- 
cal, statistical, and geometrical paradigms; “ANDING” the 
two classifiers produces highly effective results. In general, 
many classifiers could be “ANDED” but for the sonar data 
presented here, two have worked very well. Also the pro- 
cess of classifier selection relied strongly on the fact that the 
optimal feature selection process employed would be “tuned” 
to each classifier. With this in mind the two classifiers were 
chosen because: (a) they train very fast - this fact makes the 
algorithm used in the stepwise optimal selection process fea- 
sible, (b) they each divide feature space into different geo- 
metrical models — the KNN into hyperspheres and the 
ODFC into hyperplanes — this tends to make the classifiers 
complementary,(c) the classifiers use different “noise” tol- 
erant paradigms. The KNN uses sample Bayesian condi- 
tional probabilities as determined from the training data. The 
ODFC uses noise rejection and signal enhancement schemes 
based on concepts of matched-filter design. 

The KNN is a probabilistic-based neural network that 
employs radial-basis neurons [10]. The neural network has 
two layers: the attractor layer and the confidence layer. The 
attractor layer is constructed during training in the follow- 
ing way. Feature space is partitioned by a set of attractors. 
A feature vector, f, is said to belong to attractor i if 
Nf. f) <R_ where N(f. f ) is some distance measure between 
vectors f and f (for this study an L, norm was used), f is the 
center of attractor i and R_ is the radius of attractor i. 

For a given N(.,.), R., and f the probability that attractor 
i contains a given class can be estimated from sample prob- 
abilities computed from the training set. Each attractor is 
represented as a radial-basis neuron that has the feature vec- 
tor as an input and multiple outputs (one output for each 
class). The i-th neuron fires when the input feature vector 
belongs to attractor i; when it fires, the level of each output, 
associated with a given class, equals the estimated probabil- 
ity that the attractor contains that class. 

During training, the number of attractors, their center 
locations (f), radii (R,), and degree of overlap are determined 
in an optimal fashion to encompass the entire training set. 
The radii are selected using the well-known k-nearest neigh- 
bor scheme to achieve accurate probability estimates and 
adequate class boundary resolution for good class discrimi- 
nation. The more training samples in an attractor, the more 
accurate the sample probability estimates will be. In a rela- 
tive sense, an attractor with a large radius will encompass a 
large number of training samples. However, the size of the 
radii is balanced against the resolution into which feature 
space is partitioned; the smaller the radii, the better the dis- 
criminatory resolution at complex class boundaries. Finally, 
attractor overlap is used to regulate the smoothness of prob- 
ability estimates across class boundaries. 
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Training of this neural network also takes into account 
that the training data is usually statistically unfairly sampled; 
i.e., the number of training samples in a class (or subclass) 
does not reflect the true probability of occurrence of that 
class but instead reflects the artificial conditions over which 
the data was collected. For example, one might have more 
of one class than another, or more data collected under one 
condition than another. In reality, fairly sampled data is sel- 
dom available because the cost to collect an adequate data 
base is unacceptably high or simply because one does not 
have adequate access to another country’s targets or envi- 
ronment in which to collect representative data. Under such 
uncertainty, it was determined that the best statistical prop- 
erty to reinstate into the training set is that of equally likely 
classes (or subclasses). This is accomplished by a specially 
devised counting method that is used to determine the sample 
probabilities. Thus, the training of this neural network is not 
biased to favor those classes that were most plentiful in the 
training and test sets and to ignore those classes that were 
least abundant. Other neural network training algorithms 
that optimize the number of correct classifications are typi- 
cally biased by such data. 

The output of the attractor layer is then fed to the con- 
fidence layer which determines the confidence that the input 
feature vector belongs to each class. First all the attractors 
that contain the input feature vector (or are sufficiently near 
it) are identified. The confidence that the feature vector be- 
longs to a given class is then determined by interpolating the 
individual probability estimates that the feature vector be- 
longs to a class over this set of attractors. In the interpola- 
tion, the contribution of each attractor’s probability estimate 
is weighed according to how close the feature vector is to 
the attractor’s center. 

This neural network has the desirable characteristic 
that it can determine if a feature vector cannot be reliably 
classified. Specifically, if there are no attractors sufficiently 
near, then the feature vector is considered statistically dif- 
ferent than the training vectors and is said to belong to an 
unknown class. Another way to view this is to say that dur- 
ing training, the neural network never experienced anything 
that was representative of this new vector (i.e., the new vec- 
tor is not within the neural network’s experience base). Un- 
der this circumstance, there is no basis to select among the 
known classes and all outputs of the neural network are set 
to zero. Other noteworthy characteristics of this network 
are that the training algorithm is non-iterative, fast and the 
network sizes are reasonable thereby permitting very quick 
execution so specialized hardware is not needed. 

The ODFC is a classifier with its basis in linear dis- 
crimination theory [10]. Two banks of linear filters are de- 
termined from the training set: one bank is sensitive to mine 
characteristics and the other is sensitive to clutter character- 
istics. The following is a brief description on how the dis- 
crimination process works. 

Let f be a feature vector. Let F (fi) be the energy 





output for input f of the i-th filter in the mine-filter bank, and 
F (f,j) be the energy output for input f of the j-th filter in the 
clutter-filter bank. F\ and Fare linear filters; the outputs of 
these filters are linear combinations of the input features plus 
a DC offset and a warping term made up of a linear combi- 
nation of nonlinear functions of the features. The filter’s co- 
efficients, DC offset, and warping coefficients are determined 
such that, on the average, the energy output of the mine-filter 
bank is greater than the energy output of the clutter-filter 
bank when the input feature vector corresponds to a mine. 
And vice versa for a feature vector that corresponds to clut- 
ter. That is, on the average, 


max F (f(mine),i) > max F (f(mine),j) 
i j 


max F (ficlutter),i) < max F (f(clutter),j) 
i j 


The solution involves the solving of a generalized eigen- 
value problem. In the ODCF linear combinations of features 
permit boundaries between classes to be described by sets of 
hyperplanes. By adding a linear combination of nonlinear 
functions of the features one changes the hyperplane surface 
to one that is curved (warped). These nonlinear terms per- 
mit the partitioning surfaces in feature space to more effi- 
ciently fit irregular class boundaries. 

Classification is determined by “ANDING” the out- 
put of both classifiers. Both the KNN and ODFC have two 
outputs corresponding to the confidence that the input fea- 
ture vector is associated with a mine or clutter. Let C, (mine) 
and C, (clutter) designate these two confident levels for the 


knn 


KNN; and Ci a( mine) and C.4,(clutter) for the ODFC. If 


Tyan Cinn( mine) > C,,,(clutter) 
AND 
T C.,ye( mine) > C,y( Clutter) 


a mine classification is declared. The classification gain 
thresholds, 7,,, and T,,,, are adjusted to select the desired 
balance between PdPc and false alarms. 


3. Kesults 


3.1 SSS2 Example 1 


In Figures 3 through 10, Image 0 and Image 51 from 
the 60 image SSS2 data base are used to illustrate the D/C 
algorithm. Figure 3 shows the normalized image and the 
results of the detection stage. The gray scale indicator on 
the extreme right of the image indicates pixel intensity. White 
indicates a pixel intensity of 4.0 and represents the maxi- 
mum pixel intensity value displayed. Black indicates a pixel 





Figure 3 
Detections for SSS2 Image 0. 








intensity of 0.0 and represents the minimum pixel intensity 
value displayed. To give a size perspective, each sonar im- 
age has a two-meter by two-meter square highlight in the 
lower left-hand corner of the image. In all images, range 
from the sonar increases from left to right implying target 





Figure 4 
KNN Classifications for SSS2 Image 0. 
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shadows are to the right of the highlights. The detected mine- 
like targets typically resemble the target-signature mask—a 
highlight region followed by a shadow and preceded by a 
uniform background. With a strong highlight and pronounced 
shadow, Detection B is easily visible. Although Detection C 





Figure 5 
ODFC Classifications for SSS2 Image 0. 





DDFC Classification B 
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has a noticeable highlight followed by a shadow, the shadow’s 
proximity to a trough made by a fisherman’s drag net ob- 
scures the shadow region. False Detection A is triggered by 
a subtle highlight and shadow. 

Figures 4 and 5 show that both the KNN and the ODFC 
classify Detections B and C as mine-like; they both reject 
Detection A. “ANDING” the two classifiers results in Final 
Classifications B and C as illustrated in Figure 6. The squares 
in Figure 6 indicate actual mines; there one sees that both 
mines are detected and classified and one false target is de- 
tected but not classified as mine-like by either classifier. 





Figure 6 
KNN and ODFC Classifications for SSS2 Image 0. 








3.2 SSS2 Example 2 


Figure 7 shows another SSS2 example. This particu- 
lar image illustrates some of the difficulties associated with 
normalizing the background in the sonar image. Note that 
the background level of the image is not uniform. And there 
is a large band of bright pixels that extends from the upper 
left to the lower right of the image caused by a highly reflec- 
tive bottom type. Despite the lack of uniformity in the back- 
ground intensity, the AMDAC detects four candidate mines. 

In Figure 8, the KNN classifies only Detections D and 
E as mines and rejects Detections C and F. Figure 9 shows 





Figure 8 
KNN Classifications for SSS2 Image 51. 
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that the ODFC classifies Detections E and F as mines and 
rejects Detections C and D. Results from “ANDING” the 
two classifiers are shown in Figure 10. As the ground truth 
square indicates, the final mine classification is correct and 





Figure 9 
ODFC Classifications for SSS2 image 51. 
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all false alarms have been eliminated. This example clearly 
illustrates the power of combining two different classifica- 
tion approaches; separately, each approach would have had 
one correct classification and one false alarm. “ANDING” 
the two classifiers eliminated the false alarms. 





Figure 7 
Detections for SSS2 Image 51. 


Figure 10 
KNN and ODFC Classifications for SSS2 Image 51. 
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3.3 SSS2 False Alarms, Missed 
Detections and Overall 
Performance 


Overall, the AMDAC detected and classified 91% of 
the mines in the SSS2 data base with an average false alarm 
rate of 0.12 false alarms per image. Figure 11 shows the 
seven false alarms detected and three missed mines in the 
entire SSS2 data base. The false targets are quite mine-like. 
The top and bottom missed mines in Figure 11 have weak 
highlights and very poor shadows. The middle missed mine 
is embedded in very bright surface reverberation which made 
the mine signature blend in with the background. 





Figure 11 
False Alarms and Missed Detections for SSS2. 








3.4 SAS Example 


In Figures 12 through 15, Image 116 of the 215 image 
SAS image data base is used to further illustrate the D/C 
algorithm. The sonar image contains data from both the star- 
board and port sides of the sonar. Returns from the water 
column in the middle of the sonar image were blanked out 
prior to processing. Sonar operators ignore this region be- 
cause it consists of returns from the water volume prior to 
the first bottom return. Note that there are no shadows on 





Figure 12 
Detections for SAS Image 116. 














the targets; because the data was collected during shallow 
water operations, surface reverberation came into play and 
filled in the shadows. The lack of a target shadow makes the 
detection and classification process quite difficult. When a 
target exhibits only a highlight, the detection thresholds are 
typically set to accept a higher false detection rate. The clas- 
sification phase is then relied upon to reject the extra false 
targets that are detected. 

Figure 12 shows that the AMDAC has detected eight 
mine-like regions in the SAS image. All eight candidate 
regions do exhibit the desired mine-size highlight surrounded 
by a uniform background. Figure 13 shows that the KNN 
rejects two of the detections and Figure 14 shows the ODFC 
rejects five of the detections. “ANDING” the classifiers 
results in the classification of the two mines and one false 
target as shown in Figure 15. 


3.5 Kesults Summary 


Table 4 summarizes the AMDAC’s performance for 
the all three sonars, the SSS1, the SAS and the SSS2. Table 
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Figure 13 
KNN Classifications for SAS Image 116. 


Figure 14 
ODFC Classifications for SAS Image 116. 

















4 also lists expert sonar operator’s performance for the SSS1 
and SAS data bases. Operator performance for the SSS2 is 
unavailable. Because of the similarities between the SSS1 
and the SSS2, the operator performance for the SSS2 should 
be comparable to that of the SSS1. 

For the SSS1, the AMDAC detects and classifies 90% 
of the mines with 0.42 false alarms per image while the ex- 
pert operator detected 80% with 0.72 false alarms per im- 
age. For the SAS, the AMDAC performs similarly to the 
expert sonar operator at 68% PdPc. However, the AMDAC 
can also be adjusted to 92% PdPc with 0.64 false alarms per 
image. And for the SSS2, the AMDAC shows an impres- 
sive 91% PdPc with only 0.12 false alarms per image. 





Table 4. AMDAC Performance Summary 
Operator 








Fa/Image 
0.72 
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4. Conclusions and 
Future Kesearch 


Five fundamental principles have come out of this re- 
search on automated detection and classification. 

(1) The detection stage is by far the most 
computationally intensive stage. Because the entire image 
must be scanned, the basic design focus is on keeping the 
number of computations per pixel small. The main goal is 
not to detect mine-like objects, but rather to eliminate the 
majority of the image that is non-mine-like. The resulting 
detections will be processed in depth by the classification 
stages which require many more computations per pixel in 
order to extract robust and discriminatory mine-like features. 
If the detector produces 10 mine-like detections per image 
on the average (a very generous number), we have found 
that the computational requirements of our classification 
Stages are negligible compared to those of the detector. It 
then follows that this classifier’s main goal is to eliminate 
the false targets while preserving mine detections. 

(2) The classification stage uses two classifiers: the 








Figure 15 
KNN and ODFC Classifications for SAS Image 116. 














KNN and ODFC. The stepwise optimal selection proce- 
dure, described herein, dramatically improves the perfor- 
mance of each classifier by permitting a computationally 
efficient means to select the “best” subset of classification 
features from a much larger set of candidate features. As 
such it mitigates problems associated with Bellman’s curse 
of dimensionality. 

(3) Significant improvement in classifier performance 
results if the feature selection process is highly “tuned” to 
the respective classifier. Classifiers use statistical, mathemati- 
cal, and geometrical constructs to partition feature space into 
class regions. In this research we found that a subset of fea- 
tures optimal for one classifier’s partitioning scheme will be 
far from optimal for another classifier which uses a different 
partitioning approach. Therefore, we us the same partition- 
ing approach when determining the best feature subset for a 
particular classifier. 

(4) Because training procedure for both the KNN and 
ODFC is very fast, it was possible to develop a stepwise 
optimal selection procedure capable of efficiently evaluat- 
ing thousands of feature combinations. Without such an au- 
tomated procedure it was impossible for us to determine the 


“best” features. Being able to use the “best” features dra- 
matically improved the performance of our classifiers. This 
strongly suggests the following. Rapid-trainable classifiers, 
whose “best” features can be determined by stepwise opti- 
mal selection, may significantly outperform more complex 
classifiers that have computationally intense training proce- 
dures and must rely on selection methods that are either heu- 
ristic or not well “tuned” to the classifier. 

(5) There is a significant improvement in overall clas- 
sification performance that results from “ANDING” classi- 
fiers which use complementary statistical, mathematical, and 
geometrical constructs to describe class boundaries in fea- 
ture space. “ANDING” the KNN and ODFC proved very 
adept at reducing false alarms while maintaining a high prob- 
ability of mine detection and classification. This concept 
can be illustrated using the paradigm of “getting a second 
doctor’s opinion”. Two doctors look at data from a patient 
to arrive at their diagnoses. Even though they look at the 
same data, each doctor emphasizes or de-emphasizes differ- 
ent aspects based on their training and experience. If they 
conclude the same diagnosis, there is little doubt of its valid- 
ity. This also suggests that “ANDING” two or more simple 
classifiers may lead to improved performance over a single 
more complex classifier. This leads the authors to believe 
that extremely high performance in automated mine D/C can 
be achieved by “ANDING” (fusing) the diverse D/C algo- 
rithms currently being developed through ONR under their 
Mine Countermeasures research programs. 

Overall, the new detection and classification algorithm 
developed has been shown to perform significantly better 
than the detection and classification algorithm previously 
reported in SPIE 95 [1]. The new algorithm used in this 
study performs as well or slightly better than an expert sonar 
operator on the three sonar image data bases. 

The algorithm’s robustness will be further evaluated 
on new advanced synthetic aperture and side-looking sonars 
that are currently undergoing their initial sea testing. We will 
investigate the benefits of using three classifiers rather than 
just two. Using more than three classifiers produce some 
combinatorial design issues. Also, algorithms that adapt 
detection and classification thresholds will be developed in 
order to increase robustness to varying environments and 
bottom types. 
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l. IMmtroduction 


The Electro-Optic Identification (EOID) Sensors 
project is developing a Laser Visual Identification Sensor 
(LVIS) for identification of proud, partially buried, and 
moored mines in shallow water/very shallow water (SW/ 
VSW). LVIS will be deployed in small diameter underwater 
vehicles, including unmanned underwater vehicles (UUVs). 
Since the mission is mine identification, LVIS must: (a) de- 
liver high quality images in turbid coastal waters, while (b) 
being compatible with the size and power constraints im- 
posed by the intended deployment platforms. This project is 
sponsored by the Office of Naval Research, as a part of the 
AOA Mine Reconnaissance/Hunter program. 

High quality images which retain target detail and con- 
trast are required for mine identification. LVIS will be de- 
signed to produce images of minelike contacts (MLC) of 
sufficient quality to allow identification while operating in 
turbid coastal waters from a small diameter UUV. Technol- 
ogy goals for the first generation LVIS are: (a) identification 
range up to 40 feet for proud, partially buried, and moored 
MLCs under coastal water conditions; (b) day/night opera- 
tion from a UUV operating at speeds up to 4 knots; (c) power 
consumption less than 500 watts, with 275 watts being typi- 
cal; and (d) packaged within a 32-inch long portion of a 21- 
inch diameter vehicle section. 


2. Kole of EO in Mine 
Countermeasures 


Analysis of Fleet missions, including establishment of 
a Fleet Operating Area (FOA) and establishment of an Am- 
phibious Objective Area (AOA) provide a framework to de- 
termine the benefit of performing identification of MLCs 
during reconnaissance. Identification plays a key role in re- 
connaissance to establish an AOA. Reconnaissance opera- 
tions are conducted prior to the assault to help define the 
threat and to determine where mine sweeping or mine hunt- 
ing must be performed. Through reconnaissance, assessments 
are made of the extent of the mined areas, ard mine types 
and mine densities are determined. Mine identification is 
required to make these assessments. Without identification, 
only densities of MLCs may be obtained, and mine sweep- 
ing or mine hunting will be required over all portions of the 
AOA which have MLCs. With identification, mine hunting 
or mine sweeping may be restricted to the areas of the AOA 
which are mined. In typical scenarios, this can lead to a drastic 
reduction in mine sweeping/hunting assets required to com- 
plete the mission within the allotted time. 

Identification also plays a key role in reconnaissance 
prior to establishment of an FOA. Prior to moving a battle 
group or task force into an FOA, the officer in tactical com- 
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mand (OTC) must be provided an assessment of the poten- 
tial mine threat in the area. The tactical mine warfare officer 
uses data gathered by all available intelligence sources to 
estimate the likelihood and level of mining. The operational 
concept for a reconnaissance system is to use the system to 
verify the presence of a minefield in a specific area by locat- 
ing individual mines. The system is used to sample the area 
of interest, with the level of sampling effort based upon a 
balancing of the threat estimate and tactical time constraints. 
If no mines are found, the OTC either decides to accept the 
risk and continue operations or he continues sampling to in- 
crease the confidence that no mines are present. If one or 
more mines are located and identified, the battle group or 
task force will modify their intended movement for the pur- 
pose of self-protection. This will entail continuing recon- 
naissance to provide the OTC with enough information to 
decide the best course of action; i.e., either (a) map all ob- 
jects for avoidance, (b) avoid the entire minefield, (c) re- 
quest and wait for clearance by mine warfare forces, or (d) 
reconnoiter an alternative operating area. 

If reconnaissance is performed without identification, 
every MLC in the FOA must be avoided. The percent of 
operating area denied is a non-linear function of the density 
of MLCs and the radius of the avoidance circle around each 
MLC. Even quite modest densities of MLCs lead to denied 
areas which unacceptably restrict Fleet operating options. 

In summary, mine identification is an important part 
of mine reconnaissance. MLCs limit the performance of mine 
reconnaissance systems. Without identification, MLCs must 
be avoided by Fleet assets (leading to denied areas) or neu- 
tralized. Denied areas can dramatically and unacceptably 
restrict Fleet operations options, while neutralization requires 
time and appropriate Fleet assets. With LVIS providing a 
mine identification capability to a mine reconnaissance sys- 
tem, only the contacts identified as mines need to be avoided 
or neutralized. 


3. 


Impact on EO System 
Performance 


Environmental 


Environmentally derived noise sources are principal 
factors limiting the performance of underwater EO imaging 
systems. In deep, clear ocean water, an EO imaging sensor 
may have a range in excess of 100 feet, while in turbid lit- 
toral waters the same sensor’s range may be reduced to 5 
feet or less. For underwater imaging sensors, the environ- 
mental limitations may be classified as: (a) backscatter noise, 
(b) blur/glow/forward scatter noise, and (c) attenuation. Back- 
scatter noise refers to the photons which are scattered back 
into the receiver before reaching the target plane. Blur, glow, 
or forward scatter noise refers to the forward scattering of 
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photons which have been reflected from the target plane. 
Because they are scattered, they image as though they were 
reflected from the wrong position on the target, leading to a 
loss of image resolution and contrast. Finally, attenuation 
refers to the loss of photon signal due to photon absorption 
or scattering out of the field of view of the receiver. Attenu- 
ation limits signal strength at the receiver. 

A major goal of EO system designers is to optimize 
the system design to minimize the detrimental effects of the 
environment. In addition to the system architecture and pa- 
rameters, the performance of a given underwater EO mine 
countermeasures (MCM) system in a given location is de- 
termined by: (1) the Inherent Optical Properties (IOPs) of 
the water, and (2) the IOPs of the bottom and target. 

The IOPs of the water which affect the performance 
of EO MCM systems currently under development include: 
(1) the beam attenuation coefficient ¢(A)[m'], (2) the beam 
absorption coefficient a(A)[m™], and (3) the volume scatter- 
ing function B(A,O)[m"'sr'], and the scattering coefficient 
b(A)[m"]. 

The beam attenuation coefficient c(A)=a(A)+b(A) mea- 
sures the exponential decay of a narrow laser beam due to 
either absorption or scattering. The single scattering albedo 
A(A)=b(A)/e(A) measures the portion of the beam attenua- 
tion which is due to scattering. The volume scattering func- 
tion, B(A,O) measures the angular distribution of the photon 
scattering. 

The beam attenuation coefficient ¢c is the primary in- 
dicator of the distance at which an EO imaging system will 
be able to deliver acceptable imagery. Frequently, the water 
clarity is characterized according to the range R=I/c at which 
cR=1. Accordingly, when I/c=2m, for example, the water is 
said to be “2-meter water,” and when I/c=4m, the water is 
said to be “4-meter water.” Predictably, an imaging system 
will deliver acceptable imagery at a longer range — approxi- 
mately twice as far — in “4-meter water” than it will in “2- 
meter water.” The range R of an imaging system is frequently 
specified in terms of dimensionless beam attenuation lengths, 
where the range in beam attenuation lengths is given by the 
product eR. Thus, when an image is taken at a range where 
cR=5, the image is said to be taken at 5 (beam) attenuation 
lengths (AL). 

The single scattering albedo A also plays a key role in 
the performance of EO imaging systems. In typical coastal 
waters A is 70-80 percent, although it can range from less 
than 50 percent to greater than 90 percent. When the albedo 
is low, attenuation by absorption is more predominant, lead- 
ing to lower power levels. In this case, the imaging system is 
more likely to become power limited. In addition, blur/glow/ 
forward scatter is less important when the albedo is low, so 
the image contrast is higher. On the other hand, when the 
albedo is high, blur/glow/forward scatter is more important, 
leading to a reduction in image resolution and contrast. In 
addition, attenuation due to absorption is less important, so 
the system performance is less likely to be power limited. 





The volume scattering function (VSF) characterizes 
the angular dependence of the scattering, ranging from 0 
degrees (forward) to 180 degrees (backward). The VSF af- 
fects all three of the primary mechanisms of image degrada- 
tion. Backscatter depends on the magnitude of B(A,O) near 
©=n. Blur/glow/forward scatter depends on the magnitude 
of B(A,©) near O=0. Scattering at other angles contributes to 
attenuation. 

Conventional underwater imaging systems, consisting 
of a video camera and floodlights, are limited to ranges of 1- 
1.5-beam attenuation lengths (the beam attenuation length 
in coastal waters is typically 1-2.5 meters) by backscatter 
noise. They lack the range and resolution required to facili- 
tate rapid target identification. 

The EOID Sensors project has evaluated laser-based 
imaging approaches capable of reducing the impact of envi- 
ronmental noise sources in order to facilitate rapid identifi- 
cation of MLCs in turbid littoral waters at useful ranges. Two 
approaches, laser line scan (LLS) and laser range-gated 
(LRG), have been evaluated in detail with the objective of 
extending the mine identification range to five or more beam 
attenuation lengths. LLS'?’ systems reduce the effects of 
backscatter noise and blur/glow/forward scatter noise by 
synchronously scanning a narrow laser beam and a narrow 
field-of-view receiver across the sea bottom. Performance 
of LLS systems is dramatically superior to the performance 
of conventional imaging systems. LRG* imaging systems 
(i.e., imaging LIDARs) perform significantly better than 
conventional systems because they reduce the impact of back- 
scatter by temporarily “gating-out” much of the backscatter. 


4. Testing of Laser- 
Based Underwater EO 
Imaging Systems 


A test program was established in order to assess the 
relative merits of the LRG and LLS imaging approaches for 
SW/VSW mine identification. The cornerstone of this pro- 
gram was a sea test which provided simultaneous imagery 
with state-of-the-art LLS and LRG systems. This test pro- 
vided the first direct, head-to-head, comparison of the per- 
formance of these two competing technologies. This test 
endeavored to remove all questions of unequal test condi- 
tions, since the two systems were imaging the same targets 
from the same platform through the same water at the same 
location at the same time. 

The LRG system tested was Sparta Laser Systems 
Laboratory SEE-RAY LRG system. The SEE-RAY system 
employed a gated (approximately 8 nsec gate width), inten- 
sified CCD camera which generates images at standard video 
rates. F/.85 receiver optics with a 29.4 mm aperture were 
employed which gave the system an instantaneous field-of- 


view of 20 by 15 degrees. A tilt gimbal allowed movement 
of the system through an arc of 130 degrees across track, 
increasing the system field of regard. Laser illumination was 
provided by a 30 Hz, Q-switched, frequency doubled 
Nd:YAG laser operating at 532 nm with output power of 
approximately 100 mJ per 6 ns pulse. The laser output power 
could not be adjusted by the operator. However, zoom-op- 
tics on the laser allowed the output energy to be either dis- 
tributed over the receiver field-of-view or concentrated near 
the center of the image. The gating range, camera gain, laser 
zoom, and gimbal tilt angle were controlled from an opera- 
tor console. The raw LRG image data was displayed on a 
television monitor and simultaneously recorded on a 10 
MHZ, broadcast quality, U-Matic video recorder. After the 
test, digital snapshot imagery was generated from the video 
tape using an 8-bit digital frame grabbing board. Typically, 
bursts of approximately 20 sequential images were grabbed 
and visually inspected. The image deemed best was saved to 
hard disk for analysis. SEE-RAY is 39 inches long, 13 inches 
in diameter, and requires 500 watts of power. 

The LLS system tested was the ART/SAIC LS-2048. 
The LS-2048 employed a CW argon ion laser with simulta- 
neous output at 488 and 514.5 nm. Laser output power could 
be continuously adjusted from a maximum of 2 watts down 
to 65 mW in all laser lines. A rotating, four-faceted mirror 
and output optics assembly scanned the laser beam over a 
70-degree section of the sea bottom while a synchronously 
rotating four-faceted mirror and input optics assembly fo- 
cused the reflected light onto a photomultiplier tube receiver. 
The receiver circuitry incorporated a gated integrator which 
divided each 70-degree scan line into 2048 or 1024 pixels 
with an effective dynamic range of 12 bits. Raw, unproc- 
essed, 12-bit data snapshot images of 1024 x 2048 pixels or 
2048 x 1024 pixels were stored on hard disk for later analy- 
sis. Simultaneously, the scanned lines were processed as ac- 
quired and presented in a gray-scale waterfall display on the 
operator console. The continuous waterfall display was also 
recorded on S-VHS video tape. The LS-2048 is 80 inches 
long, 11 inches in diameter, and requires 5,000 watts of power. 

In order to operate in SW during daylight hours, the 
LS-2048 was modified through the insertion of 3-nm pass 
band optical filters, centered at 51 4.5 nm, into the receiver 
optics. This improvement resulted in the first successful SW 
daylight imagery using an LLS. In addition, the filter re- 
jected the LRG system's 532-nm laser light, allowing the 
LRG and LLS systems to simultaneously image the same 
targets. Of course, the 514.5-nm filter also rejected reflected 
light from the 488 line of the laser, resulting in an effective 
output laser power range (in the 514.5 line) from 800 mW 
down to 12 mW. The LLS laser powers reported in this docu- 
ment are the powers transmitted in the 514.5 line only. 

The test platform chosen for this test was Applied Re- 
mote Technology's XP-21 autonomous undersea vehicle 
(AUV) operating in a tethered mode. A payload section and 
tether were fabricated which allowed the LRG and LLS sys- 
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tems to be simultaneously installed on the XP-21. The tether 
provided power to the imaging systems; command and con- 
trol, and trickle power to the XP-21; and command and con- 
trol and image pathways to and from the imaging systems. 
Use of the XP-21 allowed precise positioning of the imag- 
ing systems, repeatable test tracks, and very efficient data 
collection. The LRG systems performance was optimized 
through allowing its operating parameters to be adjusted to 
give the best image while the XP-21 hovered above each 
target. Data runs yielding simultaneous LRG and LLS im- 
agery were typically executed only after the LRG param- 
eters were optimized by hovering over each target. 

A shallow, turbid water environment was selected for 
the test. Periodic measurements of inherent and apparent 
optical properties of the water column were made through- 
out each test day. Parameters measured included the beam 
attenuation coefficient c, the absorption coefficient a, the 
diffuse attenuation coefficient k, and downwelling and up- 
welling irradiance 

Targets deployed included positive and negative con- 
trast Air Force Resolution Targets. These targets were painted 
on 36-inch square aluminum panels using flat white and black 
paints with nominal reflectivities of 87 and 2 percent, re- 
spectively. The target designated as the “positive” target con- 
sisted of a white Air Force resolution pattern painted on a 
black background, while that designated as the “negative” 
target consisted of a black Air Force resolution pattern painted 
on a white background. This positive/negative target combi- 
nation was selected to differentiate the effect of blur/glow/ 
forward scatter from those of backscatter and attenuation by 
means of the following observations. First, notice that the 
inherent contrasts C,=( White-Black)+White of the positive 
and negative targets are identical, since they were prepared 
using the same paints on the same substrate material.* Sec- 
ond, since attenuation reduces the observed signal from the 
positive (white) and negative (black) areas by the same fac- 
tor A, it is seen that the observed contrast Cogsm(A *White)- 
(A *Black))+{A *White)=C, , cannot be reduced from the in- 
herent contrast by attenuation.° Third, since the same back- 
scatter (BS) intensity is added to the white and the black 
areas, the observed contrasts C,,.= [(A*White+BS) - 
(A *Black+BS)] + (A*White+BS) = (White-Black) + 
(White+BS+A)<C,, are decreased from the inherent contrasts 
by backscatter. However, the contrast reduction due to back- 
scatter (and attenuation) is identical for the positive and nega- 
tive targets. Fourth, white areas lose intensity to black areas 
due to forward scattering, resulting in contrast reduction. 
Since negative targets have larger white areas, they are more 
susceptible to blur/glow/forward scatter noise, and therefore 
suffer greater contrast reductions. Accordingly, greater con- 
trast reductions in negative targets are an indicator of the 
significance and magnitude of blur/glow/forward scatter ef- 
fects on image quality. 

Figures | and 2 present typical LLS and LRG images 
of the positive and negative Air Force Resolution Targets, 
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Figure 1 
Positive (white) ranges for LLS versus LRG. LLS laser power 
was 190 mW. 
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respectively. The positive (white pattern on a black back- 
ground) images in Figure | compare LLS pictures in the top 
row’ and with LRG images in the bottom row for compa- 
rable ranges as measured in beam attenuation lengths (ALs). 
The negative (black pattern on white background) images 
are similarly arrayed in Figure 2. Compare from top to bot- 
tom for LLS versus LRG contrast. These images show: (1) 
the LLS imagery is markedly superior to the LRG imagery 
for both the positive and negative targets; and (2) images of 
the bright background, negative targets are significantly 
poorer than images for the positive targets for both systems, 
but especially for the LRG system. 

The superiority of the LLS imagery is not due to any 
purported power advantage. Laser power parity for the two 
systems, in terms of the same number of photons per pixel, 





Figure 2 
Negative (black) ranges for LLS versus LRG. LLS laser power 
was 190 mW. 





LLS Range = 5.6 AL LING Range = 6.8 AL 


LAG Range = 4.7 AL LAG Range = 5.6 AL LAG Range = 6.4 AL 











Figure 3 
LLS imagery at 5.6 ALs at reduced laser powers. 
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occurred when the LLS power was approximately 200 mW. 
The LLS imagery in Figures | and 2 was taken at 190 mW 
laser power in the 514.5 nrn laser line, yielding approximate 
energy parity with the LRG system . Figure 3 shows compa- 
rable LLS imagery taken at 5.6 ALs using 95, 35, and 12 
mW laser power in the 514.5 laser line -- far below the en- 
ergy density employed by the LRG system. Quantitative 
analysis reveals that the LLS contrasts are maintained at these 
low power levels, even though shot noise increases.’ Con- 
versely, when the LRG system laser was "zoomed down" on 
the targets, resulting in an approximate four-fold increase in 
energy density on the target, there was no observable im- 
provement in LRG image quality. 





Figure 4 
Contrast versus range plot for LLS and LRG images. 
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Contrast values were extracted from various three-bar 
patterns on these targets as a quantitative assessment of the 
relative quality of the images. Figure 4 presents these con- 
trasts as extracted from the vertical three-bar pattern above 
the large square (e.g., the bar number 2 of series 0). The 
widths of these bars are 2.24 cm, while the sides of the large 
square are 10 cm. The sea test data has been supplemented 
with data from an earlier LRG tank test. Figure 4 shows the 
contrast points for the LRG and LLS systems are quite dis- 
tinct and well separated. At any given range, the observed 
contrast of the LLS system is typically two to three times the 
contrast of the LRG system. Moreover, in order to match 
LLS contrasts, the LRG system frequently needed to be two 
or three times closer to the target. 

The imagery and the contrast graph clearly show that 
the contrasts of the negative targets are significantly lower 
than the contrasts of the corresponding positive targets. As 
indicated above, this disparity between the positive and nega- 
tive target contrasts cannot be explained by means of attenu- 
ation and backscatter. Rather, it is a direct indicator of the 
dominant role of blur/glow/ forward scatter noise in the deg- 
radation of the image quality. This dominant effect is par- 
ticularly evident in the LRG system, where the negative tar- 
get imagery is far inferior to the corresponding positive tar- 
get imagery. For the LLS system, the negative target con- 
trasts are smaller than, but still reasonably close to, the cor- 
responding positive target contrasts, indicating success for 
the LLS system in partially rejecting blur/glow/forward scat- 
ter noise. 

Even more direct evidence of the dominant role of blur/ 
glow/forward scatter noise in the degradation of LRG imag- 
ery exists. Prior to the sea test, a test of a breadboard LRG 
system had been completed using a 4- x 4- x 90-foot test 
tank. The breadboard LRG was a dry unit; a large periscope 
was used to insert the laser pulse into the tank and to extract 
the return image pulse from the tank. This arrangement al- 
lowed data collection in a novel, "back-illuminated" con- 
figuration in addition to the normal "front-illuminated" con- 
figuration. in the back-illuminated configuration, backscat- 
ter noise was completely non-existent, allowing a very 
straightforward evaluation of the role of blur/glow/forward 
scatter noise in LRG systems. Figure 5 is a diagram of the 
tank test configuration showing the two laser paths. This test 
was conducted using fresh water, with Maalox used to add 
controlled amounts of turbidity to the water. Tests were con- 
ducted under simulated turbid, littoral water conditions, with 
beam attenuation lengths of 3, 2, 1, and 0.5 meters. 

In the front-illuminated configuration (path 1), the la- 
ser beam is inserted into the water using a periscope and 
propagates through the water on its way to the target. The 
image pulse is reflected from the target, propagates through 
the water on its way back to the periscope, where it is ex- 
tracted from the water and recorded on the gated camera. In 
this configuration, the image is degraded by all of the envi- 
ronmental noise sources; i.e., attenuation, backscatter and 
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Figure 5 
Test tank facility. 
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power-dependence of the image. Targets used in the back- 
illuminated configuration mode included positive and nega- 
tive Air Force resolution patterns silk-screened on transpar- 
ent plexiglass sheets, and a "point source" target. The point- 
source target consisted of six holes (diameters 8.33, 6.4, 4.0, 
2.8, 2.0, and 1.0 mm) drilled in two rows (the two largest 
and the four smallest) in a sheet metal plate. The rows were 
separated by 232 mm, while the holes were separated by 51 
mm. 

Figures 6 and 7 show typical data sets from this tank 
test. Figure 6 illustrates LRG tank test data at a range of 3- 
beam attenuation lengths with front and rear illumination of 





Figure 7 
LRG tank test images 5.6-beam ALs, normal and high power. 





blur/glow/forward scatter noise. Contrast data from this con- 
figuration has been included in Figure 4. 

In the back-illuminated configuration (path 2), the la- 
ser pulse propagates over the tank to the back side of the 
target, is inserted into the water by a second periscope, passes 
through a diffusion apparatus before back-illuminating trans- 
parent targets. In this configuration, backscatter noise is com- 
pletely non-existent. In addition, the laser beam does not 
suffer attenuation due to scattering and absorption by water 
on its way to the target, so the intensity of the image leaving 
the target plane can be considerably higher than the corre- 
sponding image in the front-illuminated configuration. Ac- 
cordingly, much of the back-illuminated imagery was col- 
lected in two modes. In the first mode, the laser pulse was 
attenuated so that the laser power leaving the target was ap- 
proximately the same as the corresponding laser power in 
the front-illuminated configuration. In the second mode, the 
laser pulse was not attenuated, allowing an evaluation of 





Figure 6 
LRG tank test images, 3-beam ALs. 
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the positive and negative targets. The outline of the 30-inch 
diameter “back-illumination spot" is clearly evident in the 
back-illuminated positive and negative target images. The 
lower right hand panel shows imagery of the back-illumi- 
nated “point-source” target. Figure 7 provides additional LRG 
tank test images at 5.6 ALs. The two top panels are standard 
front-illuminated, positive and negative targets. The panels 
in the middle line are back-illuminated, positive negative, 
and point targets, all at normal laser power. "Normal laser 
power” means that the back-illuminating laser power was 
attenuated to approximately match the laser power incident 
on the front-illuminated targets. (These powers were mea- 
sured with a radiometer). The Figure 7 bottom three panels 








Figure 8 
Contrast versus range for LRG tank test. 
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are the same as the middle row, except that the back-illumi- 
nating laser power was not attenuated, yielding a laser power 
leaving the target significantly higher than the "normal laser 
power" leaving the front-illuminated, reflective targets. In 
this case the laser power was approximately 25 times higher 
than normal laser power. 

Figure 8 presents contrasts extracted from this data for 
the front-illuminated, positive and negative Air Force Reso- 
lution Targets and corresponding back-illuminated positive 
and negative resolution targets for the back-illuminated tar- 
gets. The contrasts of the back-illuminated targets are, as 
expected, higher than that of the corresponding front-illumi- 
nated targets. However, the contrast of the front-illuminated 
positive targets is significantly higher than the correspond- 
ing contrasts of the back-illuminated negative targets. This 
data places significant bounds on the magnitude of the con- 
trast reduction due to backscatter; if the contrast reduction 
of the LRG system was mainly due to backscatter, the con- 
trasts of the back-illuminated negative (and positive) targets 
should be significantly higher than the contrasts of the front- 
illuminated positive (and negative) targets. Moreover, the 
contrasts of the back-illuminated positive and negative tar- 
gets should be close to identical. Neither of these is the case. 

The data from the "point source” target clearly dis- 
plays the fact that the LRG system was not power limited 
under the conditions of this test. The LRG system had enough 
sensitivity and gain to image the point source target at 5.6 
ALs; it is blur/glow/forward scatter which limits the system 
performance, particularly for the negative targets. In fact, 
the system had enough sensitivity and gain to resolve the 
8.33- and 6.4-mm point sources at 9-beam attenuation 
lengths. Under these conditions blur/glow/ forward scatter 
noise rendered even the positive back-illuminated image fea- 
tureless. 

It should also be noted that the LLS imagery and con- 
trasts from the sea test is far superior to even the back-illu- 


minated, enhanced power, LRG tank test imagery and con- 
trasts. 

Shortly after the completion of the sea test, a second, 
very brief, tank test was conducted which compared the per- 
formance of the LRG system in the sea test configuration 
with a prototype, enhanced receiver, LRG system. The pro- 
totype receiver has a higher quantum efficiency photocath- 
ode (approximately 45 percent versus approximately 15 per- 
cent) and a lower excess noise factor (approximately 1.1 
versus approximately 3). However, these system improve- 
ments did not lead to significantly improved imagery under 
turbid littoral water conditions. In fact, under most condi- 
tions, the sea-test configuration performed at a par with, or 
superior to, the prototype configuration. Because this was a 
prototype receiver, special caution must be exercised when 
drawing conclusions. However, because the system was 
tested under turbid water conditions, the system was not 
power limited. Rather, its performance was limited by blur/ 
glow/forward scatter noise, which was not addressed by these 
receiver enhancements. 

The positive and negative Air Force Resolution Tar- 
gets provide examples of imagery of high inherent contrast 
targets. Because of their high contrast, they are very coop- 
erative targets. Caution must be exercised when extrapolat- 
ing these results to imaging less cooperative, low inherent 
contrast targets. 


5. 


Conclusions 


This completed the first direct, head-to-head, perfor- 
mance comparison of LLS and LRG underwater imaging 
technologies. The results were an unambiguous appraisal of 
their relative merits for the mine idertification mission in 
turbid, littoral waters. The following conclusions were 
reached: 

The LRG system demonstrated: (a) successful detec- 
tion (but not identification) of high contrast targets at rela- 
tively long range, (b) insensitivity to ambient light, and (c) 
insensitivity to platform motion. However, it had extreme 
difficulty identifying or even detecting low contrast targets 
even at relatively short ranges. It was found that the contrast 
and resolution required for the mine identification mission 
was severely limited by blur/glow/forward scatter noise. 

The LLS system demonstrated the contrast and reso- 
lution required for mine identification under a robust set of 
conditions, including: (a) successful identification of high 
and low contrast targets, (b) successful imagery at high and 
low laser power, (c) successful imagery of tethered targets 
in the volume, and (d) successful imagery during day and 
night conditions. The system demonstrated effective reduc- 
tion of blur/glow/forward scatter noise, as well as backscat- 
ter noise. Under certain situations, the LLS system became 
power limited, suggesting improvements in detector quan- 
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tum efficiency and increased laser power would be benefi- 
cial. 

Based on these test results, requirements analysis, and 
sensor performance modeling, the decision has been made 
to base the LVIS sensor on LLS technology. The LVIS sen- 
sor has been designed, and is currently being fabricated by 
Applied Remote Technology. It is scheduled for initial sea 
testing in July/August 1995. 
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Other definitions of contrast are possible, and may be 
more useful for certain puiposes. This definition was 
used because it supports the separation of the effects 
of attenuation and backscatter from those of blur/ 
glow/forward scatter. 

Of course, if the signals are sufficiently attenuated, 
the image will be lost in the shot noise. Sufficient 
power and gain is required. However, shoi noise is 
not the primary issue for the images under consider- 
ation. 

The LLS images exhibit image distortions caused by 
non-uniform platform motions, especially roll 
induced by propeller torques when the platform is 
accelerating. 

These images were taken during early evening, when 
ambient light levels were low. Comparable images 
taken near noon, when ambient light levels were 
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high, led to significant contrast reductions at the 
lowest laser power levels due to ambient light effects. 
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Abstract 


Magnetic sensors offer a complementary approach to 
active acoustics for shallow water mine reconnaissance and 
hunting. The U.S. Navy has developed an approach using a 
5-channel tensor magnetic gradiometer to provide enhanced 
classification and location in mobile operations beyond the 
capability offered by the commonly used single channel to- 
tal field magnetometer. Buried mine detection and low false 
alarm rates were demonstrated using such a 5-channel gra- 
diometer in fusion with acoustic sensors. This sensor fea- 
tured bulk and wire niobium low critical temperature (low 
T.) superconducting components cryocooled by liquid he- 
lium. Advances in material research and new concepts are 
being pursued to enhance opportunities with this 5-channel 
gradiometer concept: an advanced low T_ superconducting 


prototype incorporating all thin film niobium superconduct- 
ing components to demonstrate increased detection range in 
the Joint Countermine Advanced Concept Technology Dem- 
onstration, a high T_ superconduting concept using liquid 
nitrogen refrigeration to reduce package size and cryogenic 
support requirements, and a room temperature fluxgate pro- 
totype for man-portable applications where shorter detec- 
tion ranges are useful. Recent experimental results using 
this technology have been obtained to demonstrate an en- 
hanced capability for the detection of unexploded ordnance 
for environmental cleanup. In this paper, these recent ad- 
vances in sensor development and the new testing results 
will be reviewed. 
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I. Introduction 


Magnetic sensors have proven merit for mobile area 
surveys and search operations conducted from air, land or 
sea including application for the detection, classification, and 
localization of sea mines, unexploded ordnance (UXO), and 
chemical, biological and nuclear waste [1]-[4]. In fact, in an 
assessment conducted by the Jet Propulsion Laboratory for 
the Army Corps of Engineers, magnetic sensors in general 
were identified among the most useful sensors for UXO de- 
tection and localization, and superconducting gradiometers 
were specifically identified as the most useful tool in a class 
by themselves [5]. For applications in sea mine countermea- 
sures, we Can envision operational scenarios in which long- 
range detection is required for reconnaissance and hunting 
in preparation for an amphibious assault or shorter-range 
detection is required for diver mine detection and avoidance 
(Fig. 1). 

In Sections I.A. and I.B., we shall describe two types 
of sensors which detect magnetic anomalies: sensors which 
detect changes in the local magnetic field, magnetometers, 
and sensors which measure the spatial derivatives of mag- 
netic field, (first order) gradiometers. 


A. Magnetometers 


The performance of a magnetic sensor is measured by 
its detection range, which is a function of its configuration 
and sensitivity and the magnetic moments of the targets of 
interest. In the far field, a target can be well approximated as 
a magnetic dipoie. In this approximation, relatively simple, 
analytic expressions can be written to relate sensitivity re- 


quirements in terms of nominal values for target magnetic 
moment and range. 

For the special case in which a circular, connected con- 
ducting loop with area A is carrying an electrical current I, 
we can define the magnetic moment of the loop as the vector 
m = IAf where fi is the normal to the loop in the direction 
defined by the right-hand rule for positive current. The In- 
ternational System (SI) unit for magnetic moment is ampere- 
meter-squared (A-m?). The magnetic induction variations in 
the Earth’s field, without the use of very sophisticated for a 
magnetic dipole can be written as 


p= fe| Amr 2) 
4n 


r° 7 

The concept of magnetic moment m can be general- 
ized for an arbitrary magnetic body and the magnetic induc- 
tion of the body will approach the result of (1) in the far field 
In air, the magnetic induction B is related to the magnetic 
field H via B=y,H. Hereafter, B will be used exclusively 
and will be referred to as the magnetic field. The following 
units for magnetic field will frequently be used to represent 
sensor sensitivities: nanotesla (1 nT=10°T), picotesla 
(1 pT=10"T), and femtotesla (1 fT=10°°T). 

A number of sensors, notably fluxgate magnetometers 
and superconducting magnetometers, measure the individual 
vector components of field. A three-axis vector magnetom- 
eter, likely using fluxgate or superconducting sensors, is very 
useful for localization (providing three channels of informa- 
tion). For stationary applications in geophysics and barrier 
defense, such sensors are effective. However, to date, these 
sensor types have not proven effective for mobile applica- 
tions, since means to compensate the anomalous signals aris- 





Figure 1 


Operationals scenarios for the magnetic detection and classification in sea mine countermeasures: (a) operation onboard 
an unmanned underwater towed vehicle towed behind a semi-submersible remotely-operated vehicle for reconnaissance 
and hunting and (b) diver-portable operation for diver mine detection and mine avoidance. 
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ing from rotations in the Earth’s magnetic field have not been 
devised. 

Other magnetometers, notably those based on nuclear 
or atomic resonance processes, measure the magnitude of 
the total magnetic field and are known as total-field magne- 
tometers. Let B, denote the magnetic field of the Earth, and 
let b denote the induction generated by an anomaly. If 
lbl << IB |, then the signal observed by a total field sensor 
(referenced to the baseline Earth field) is 





b. = |B, + b|-|B,| = /B +2B, eb+b° «-s ; 2) 
0 


As a result of the right-hand approximation above, a total- 
field magnetometer does not simply measure the magnitude 
of the magnetic-field anomaly, but measures instead the pro- 
jection of that anomalous signal onto the earth’s field 

Total-field magnetometers have provided the gener- 
ally accepted method for magnetic anomaly detection. In 
particular. the AN/ASQ-81 and its successor the AN/ASQ- 
208 are total-field magnetometers utilized by the U.S. Navy 
for submarine detection from airborne platforms. A major 
advantage of this type of sensor is its insensitivity to rota- 
tion in the Earth’s background field of 50,000 nT (since total 
field is a rotational invariant). 

Measurements by total field magnetometers are diffi- 
cult to interpret because these sensors effectively measure 
the projection of the anomalous magnetic field vector onto 
the Earth’s magnetic field instead of the total field. Interpre- 
tation often requires an experienced operator, and precise 
anomaly locations ate difficult to obtain. Since total field 
magnetometers provide only one channel of information, they 
lack valuable target vector information. In particular, they 


provide very limited localization and little capability for 
anomaly classification through moment determination. More- 
over, these sensors are limited in field operation to sensitivi- 
ties at levels approaching 0.1 nT as a result of geomagnetic 
noise, i.e., temporal variations in the Earth’s field, without 
the use of very sophisticated compensation schemes. 


B. Gradiometers 


The gradient of the magnetic field (in standard MKS 
units of T/m) is a second-order tensor with components given 
in Ref. 6 by 


OB 
G pated 5 
ij ox, 


= - Heme rn -1°6,)—r? (rm, + rm,)} (3) 


As a result of Maxwell’s equations in free space, only 5 of 
these 9 tensor elements are independent. For this reason, (first 
order) gradiometers are typically designed with 5 indepen- 
dent gradient channels, using the minimum number which 
permits characterization of the local tensor gradient field. 

It is feasible to determine the bearing vector and the 
magnetic-moment vector direction of the dipole by inver- 
sion of the gradient equations at a single point only [6]. More 
recently it has been shown that the addition of gradient rate 
information at a single point leads to a unique solution for 
dipole position and moment vector [7]. 

The contradiction of the gradient tensor defined by 





Figure 2 


Gradiometer sensor concepts including (a) a single SQUID-based gradiometer channel and (b) a conceptually simple 5- 
channel gradiometer configuration capable of magnetic dipole localization and moment classification (with 3 orthogonal 


magnetometers for motion compensation). 
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is a rotational invariant associated with the gradient tensor 
analogous to the magnitude of the field vector. This quantity 
may prove very useful for applications in which a gradiom- 
eter is subjected to large rotations during the period of mea- 
surement, e.g., hand-held operation (in contrast to straight 
runs onboard stabilized platforms) [6]. 

An example of one configuration to measure a single- 
gradient tensor component and a simple configuration to 
measure 5 independent gradient components are displayed 
in Figure 2. Each gradient tensor component is measured by 
a spatially separated loop pair connected in a common-mode 
rejection configuration. 3 vector magnetometers are included 
in the 5-channel gradiometer displayed in Figure 2(b) to com- 
pensate for the residual magnetometer signals in the gradi- 
ometer channels that arise as a result of manufacturing im- 
perfections in the gradiometer loops. 


Gradiometers offer the potential to remove many of 


the limitations associated with magnetometers because the 
output of a gradiometer is typically produced by twin mag- 
netometers operating in differential mode. In particular, this 
configuration provides common-mode rejection of the nomi- 
nal 0.1-nT temporal variations in the Earth’s field and of the 
nominal 1000-nT field changes arising from typical 1 de- 
gree sensor rotations while in towed motion. 

Gradiometers may be fabricated using many available 
magnetometer technologies. Available fluxgate and total-field 
magnetometers can perform at levels approaching 1-10 pT, 
while superconducting magnetometers utilizing Supercon- 
ducting Quantum Interference Devices (SQUIDs) can per- 
form at levels on the order of 1-10 fT. The extreme perfor- 


mance available from superconducting magnetometers pro- 
vides a capability to fabricate gradiometers with high sensi- 
tivity having short baselines. The short baseline provides a 
compact package with extreme coherence between the mag- 
netometers. The coherence is required to maintain high per- 
formance when in motion, and the compact package is ame- 
nable to implementation aboard one vehicle. In addition, with 
the compact sizes, a large number of gradient channels can 
be integrated into a small package to obtain complete posi- 
tion and moment determination of a magnetic dipole target 
at a single point in space. 


C. The Superconducting 
Gradiometer/Maqnetometer 
Sensor 


Almost all of the efforts with SQUID sensor technol- 
ogy have dealt with sensors inside a very controlled labora- 
tory environment, and to a more limited extent outside the 
laboratory at stationary locations, notably for geophysical 
measurements. During the | 980’s, the Coastal Systems Sta- 
tion (CSS) developed the Superconducting Gradiometer/ 
Magnetometer Sensor (SGMS) specifically for mobile op- 
erations outside the laboratory environment (Fig. 3). The 
sensor employs largely niobium bulk and wire superconduct- 
ing technology (with thin-film Josephson Junctions), features 
dc SQUIDs housed inside superconducting shields, and is 
convectively cooled to 4 degrees Kelvin by helium gas evapo- 
rating from a liquid helium reservoir. Under some field con- 
ditions, the SGMS has attained sensitivities on the order of 
lpT/m-Hz" at 0.1 Hz. In comparison total-field gradiom- 
eters (with a 0.3-meter baseline) specifically designed for 
operation onboard a moving platform have achieved sensi- 





Figure 3 


Major subassemblies in the SGMS package: (a) dewar and (b) sensor probe unit. 





(a) 
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tivities on the order of 30 pT/m-Hz'” at 0.1 Hz and a fluxgate 
gradiometer described in Section VII.A has achieved a sen- 
sitivity on the order of 300 pT/m-Hz"™ at 0.1 Hz [1], [2]. 
The CSS initiated the Magnetic and Acoustic Detec- 
tion of Mines (MADOM) Project starting in 1985. This 
project successfully demonstrated the value of magnetic and 
acoustic sensor fusion for mine reconnaissance and hunting. 
The SGMS demonstrated, for the first time, high sensitivity 
and rugged, robust, and reliable performance of a supercon- 
ducting gradiometer operating outside the laboratory envi- 
ronment onboard a towed underwater vehicle with sea test- 
ing conducted for a period of 7 years. Gradiometer opera- 
tion was automated (with the exception of semi-automated 
initial tuning) and fully-automated, real-time magnetic de- 
tection and classification signal processing was demonstrated 
to provide effective and accurate moment determination and 
localization for single and multi-target cases [1], [2]. 


D. Quantitative Comparison of 
Magnetometers and 
Gradiometers 


The signal strength of a magnetic dipole decreases as 
the third power of the range for magnetic fields and as the 
fourth power of the range for magnetic field gradients. The 
approximate ranges of magnetometers and gradiometers are 
displayed in Figure 4 as functions of dipole strength and sen- 
sor sensitivity. It can be shown that the sensitivity require- 
ments for a magnetometer and a gradiometer, respectively, 
to have the same detection range r against a given dipole 
target, is given by the approximate relation NIN, ~ 3/r. For 





Figure 4 

Approximate ranges (in meters) of magnetometers and 
gradiometers as a function of target strength in terms of 
magnetic moment (in units of A-m*). Curves are given for 
different sensor sensitivities (in units of nT for magnetometers 
and nT/m for gradiometers assuming a 10 dB signal-to-noise 
ratio in both cases). 
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Figure 5 

A comparison of magnetometer and gradiometer capabilities 
for target localization. The profiles measure tf = total field 
(Equation (2)) and gr = gradient magnitude (Equation (4)). 
Profiles are given for a magnetic dipole centered at origin 
with (a) a general orientation and (b) orientation in an east- 
west direction. Observe that the gradiometer profiles are 
approximately circularly symmetric about the dipole’s location 
so that “gradient” searches normal to the gradient magnitude 
profiles are meaningful for the gradiometer. In contrast, the 
magnetometer profiles are not amenable to such 
straightforward interpretation. 
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example, the detection of a mortar shell with a magnetic 
moment of 0.1 A-m/ at a range of 3 m, requires a magnetom- 
eter with sensitivity of 0.36 nT or a gradiometer with sensi- 
tivity of 0.36 nT/m (given a 10-dB signal-to-noise ratio for 
both cases). 

It should be noted that the high rate of signal reduc- 
tion with the fourth power of distance in the case of a gradi- 
ometer represents an apparent shortcoming for a gradiom- 
eter configuration. We believe that the ability to develop gra- 
diometers with sensitivity greater than 1x10° nT/m for mo- 
bile operations and the extreme difficulty in utilizing mag- 
netometers with sensitivity greater than 0.1 nT in mobile 
operations significantly outweighs this shortcoming. More- 
over, the fourth power reduction of detection range with 
moment for a gradiometer has merit for the detection of tar- 
gets with relatively small moments. For example, a 3 pT/m 
gradiometer can detect an individual 500-pound bomb (with 
a moment of 30 A-m7’) at a range of 33 m and a 60-mm mor- 
tar shell (with a moment of 0.1 A-m*) at a range of 8 m. 
Hence there is only a factor of approximately 4 reduction in 
detection range for a 60-mm mortar shell compared to the 
500-pound bomb although there is a factor of 300 in reduc- 
tion of magnetic moment. 

An example of the ease of interpretation for 5-channel 
gradiometer data compared to single-channel total-field 
magnetometer data is displayed in Figure 5. Magnetic pro- 
files have been generated for a 60-mm mortar shell buried | 
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meter under ground for two different orientations with re- 
spect to the Earth’s background field. In this example, the 
magnetometer profiles and gradiometer profiles are given 
by the anomalous total field, Eq. (2), and the corresponding 
changes in magnitude of the gradient tensor, Eq. (4), respec- 
tively. The complex total field profiles require precision data 
and critical interpretation to localize dipole sources. The sym- 
metric gradiometer profile leads to straightforward interpre- 
tation convenient for gradient searches for dipole localiza- 
uion. 


ll. Scenarios for 


Gradiometer Mobile 
Operation 


We can envision three general types of operational 
scenarios: relatively long-range rapid surveys for target clus- 
ters, more moderate-range searches against individual tar- 
gets, and detailed close-in surveys. The selection of mag- 
netic sensor type will largely depend on these operational 
requirements, determined primarily by the desired detection 
range, which is, in turn, a function of the magnetic moment 
of the targets and of sensor sensitivity. The selection also 
depends on such factors as financial budgets, logistical sup- 
port and technical expertise of the operators. To date, mag- 


netic sensor approaches have provided limited localization 
and mapping capabilities. To gain widespread acceptance, 
approaches must be introduced which provide accurate lo- 
calization and target classification, and which lend themselves 
to straightforward interpretation and minimal training. Per- 
formance must not be limited by magnetic noise from the 
host platform and other subsystems. For land-based opera- 
tions, the system must be capable of operating over rough, 
overgrown terrain. The sensor and associated signal process- 
ing also must deal effectively with environmental noise. 


A. General Considerations 


Long-range rapid surveys conducted from an aircraft 
have been proposed for initial surveys to locate clusters of 
UXO targets [4]. The mine reconnaissance/hunting demon- 
stration of MADOM represents an example of a moderate- 
range search. Land-based manual surveys or diver opera- 
tions for mine detection or avoidance provide exampies of 
close-in surveys. High sensitivity will be critical for the long- 
range and moderate-range scenarios, but such sensitivity will 
likely be sacrificed for the close-in surveys at ground level. 
In fact, local geologic noise limits a gradiometer’s noise floor 
to levels on the order of 0.05 up to 10 nT/m-Hz'” which may 
limit the use of high-sensitivity gradiometers such as the 
SGMS for some land-based operations. 

A perspective on the role of higher sensitivity gradi- 





Figure 6 


Two modes of operation for 5-gradient channel sensors: (a) wide-area surveys using high sensitivity sensors and (b) short- 


range searches for single targets using less sensitive sensors. 
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ometers used for wide-area searches and lower sensitivity 
gradiometers for close-in surveys can be obtained from the 
following example (Fig. 6). A 3 pT/m gradiometer can de- 
tect a grouping of twenty 500-pound bombs (clustered in a 
circle several meters in diameter) at a range of 46 meters. An 
area search rate of 1 km*/hr can be obtained when the sensor’s 
altitude is 15 meters moving at a forward speed of 15 km/hr. 
When deployed from the ground, a less sensitive 300 pT/m 
gradiometer would provide detection ranges of 10 meters 
and 3 meters for the detection of a 500-pound bomb and a 
60-mm mortar shell, respectively. 

Two fundamental approaches for operation of a 5-chan- 
nel tensor gradiometer stand out. First it will often be conve- 
nient to conduct straight-line searches at constant velocity 
for airborne and underwater vehicle operation and also for 
some land-based surveys if the terrain admits straight-line 
trajectories. For other cases, including many man-portable 
land-based and diver operations, it will be unreasonable to 
expect controlled operator motions along straight-line tra- 
jectories at constant velocities. In either case, we can imple- 
ment signal processing using point-by-point inversion of the 
motion-compensated gradiometer signals. For the special 
case of straight-line searches, we can also utilize least-squares 
fit to the time series in order to improve target detection. 
One approach using point-by-point inversion is described in 
Section II.B and one approach using a least-squares fit is 
described in Section II.C below. 


B. Localization using Point-by- 
Point Signal Inversion 


As mentioned above, the five independent gradient ten- 
sor components at a single point can be used to construct the 
bearing vector to a dipole source, and a scaled moment vec- 
tor with the same direction as the dipole moment vector and 
a magnitude given by m/r* [6]. The difficulty with this inver- 
sion is that there are multiple solutions: two nontrivially re- 
lated solutions in a given half space, and two additional so- 
lutions obtained by reflection of the first two through an ori- 
gin centered on the gradiometer. This multiplicity of solu- 
tions has limited the practical application of this algorithm. 
More recently, Wynn [7] has investigated the use of the rate 
of change of the gradient tensor components. and their role 
in resolving the scaling and uniqueness issues associated with 
the gradient tensor inversion. 

For a sensor with specified translational velocity v, the 
time rate of change of the gradient tensor components has 
the form 


G, = Foal 35(mer)(ver)r,r,- Sr?[(mev)rr, 


+(mer)(vr+vir,) + (ver)(mri+m r) 
+(mer)(ver)5,] + r[my,+ mv, + (mev)8,J) (5) 


These equations have been inverted to give multiple 
solutions for the bearing vector and a scaled moment vector 
in the direction of the dipole moment vector with a magni- 
tude given by m/r°. The solution for bearing vector and mo- 
ment vector direction common to the two inversions is 
unique, and the different scaling for the scaled moments in 
the two inversions yields the range to the dipole, resulting in 
a unique solution for m and r on a point-by-point basis. In 
practice, it is not necessary to specify a sensor velocity. All 
that is needed is knowledge of the position of the gradiom- 
eter relative to the Earth’s reference frame. Work is ongoing 
to apply these algorithms to the practical interactive local- 
ization of buried UXO by means of a man-portable tensor 
gradiometer. 


C. Localization using a Time- 
Series Least-Squares Fit 


A mathematical model for detection, classification, and 
localization (D/C/L) of multiple stationary magnetic dipole 
targets using a gradiometer (with 5 independent tensor gra- 
diometer channels appropriately selected and 3 orthogonal 
magnetometer channels) moving in a straight line trajectory 
past the targets at a constant speed has been developed and 
validated in the MADOM project. The 5 output signals S(t) 
(i=1,2,....5) from the 5 independent gradiometer channels and 
the 3 magnetic field components B (t) (I=1,2,3) are measured 
as a function of time as the sensor moves past the targets. 
The time derivatives of field dB(t)/dt are calculated from 
the B (t) for eddy-current compensation. The objective of 
this model is to extract the dipole signals G,((k=1,2,....m) 
for the unknown number n of dipole targets and then to de- 
termine the magnetic moments and the positions of the n 
targets. 

The model describes the signals S, in the 5 gradiom- 
eter channels by the equations 


S,(t)= Ye, {Soo +0, + ¥ B,B,(*) + p ¥,B,(t)+v,(t) 


(6) 


where c is the pre-determined calibration matrix for 
the gradiometer, v,(t) and a, are uncompensated noise (set- 
ting the noise floor of the channel compensation parameters 
per channel) and channel biases, respectively, and B,, and y, 
are the balance and eddy-current vectors for channel i. An 
iterative analysis first estimates the a’s, B’s, y's, (a total of 
35 parameters) and then executes a gradient search for the 
location of the single target that best fits the residual signal. 
The a’s, B’s, y's, and target location and moment are then 
optimized, and the procedure is repeated for a second target. 
Targets continue to be added until finally no target can be 
found whose signal contributes substantially to a reduction 
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Figure 7 
The five gradiometer signals in a 45-m section of data on a 
pass near a Mk83 1000-lb bomb. 
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in total signal power, at which point the algorithm termi- 
nates. 

Fig. 7 displays the motion-compensated signals ob- 
tained from the gradiometer and the information extracted 
from the algorithm in detecting a 1000-Ib Mk83 bomb. The 


five gradiometer signals are displayed for a 45-m section of 


data with the bomb located 15 m to port and 5 m below the 
gradiometer at the closest point of approach. The measured 
and predicted target location, moment, and orientation are 
also displayed to indicate the capability of the algorithm to 
provide good fits to the data. 


il. Mobile Underwater 
Debris Survey System 


A project has been initiated to develop and evaluate a 
Mobile Underwater Debris Survey System (MUDSS) capable 
of finding and accurately mapping the locations of UXO rang- 
ing from small shells to large bombs in water depths of from 
4 to approximately 100 feet in coastal regions at formally 
used defense sites [3]. The effort involves a collaboration 
between CSS and the Jet Propulsion Laboratory [8]. In addi- 
tion to the application of UXO detection at coastal sites uti- 
lizing underwater towed sensor suites as described in this 
section, the use of a superconducting gradiometer in a sys- 
tem concept for rapid airborne reconnaissance and survey of 
UXO sites has been proposed [4]. 


A. General Project Description 


The MUDSS project is divided into two phases. Phase 
I, which ran for one year and culminated in an at-sea feasi- 
bility demonstration of a multi-sensor MUDSS prototype 
against UXO in a drill target field The feasibility demon- 
stration was successfully executed in August and September 
of 1995 in St. Andrews Bay (near Panama City, FL) with the 
SGMS utilized for magnetic detection. Phase II will culmi- 
nate in 1997 in a technology demonstration consisting of a 
UXO survey at a yet-to-be-determined formally used defense 
site. 

Figure 8 depicts the first version of the MUDSS sys- 
tem which was fielded in 1995 for the feasibility demonstra- 
tion. The surface craft is a custom designed magnetically 
and acoustically quiet. shallow draft, trailerable catamaran. 





Figure 8 


The Mobile Underwater Debris Survey System: (a) artist's concept detailing key features of the system and (b) photograph 
of system in transit to test site with the dead-weight depressor housed out of water for speed and manueverability. 
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Figure 9 


Gradiometer targets found in one run over the linear field with the open circles ° indicating the actual target locations and the 
solid circles + indicating the positions predicted by the D/C/L algorithm. The size of the solid circles + indicates the predicted 


magnitude of the target’s magnetic moment. 












































































































































A dead-weight depressor is suspended off of the back of the 
catamaran to maintain the appropriate depth for sensor op- 
eration and to house part of the sensor suite (a RESON Seabat 
ahead looking sonar, a CSS-developed high frequency 
sidescan sonar, and a Raytheon-leased LS 4096 laser 
linescanning electro-optic sensor. A second neutrally buoy- 
ant tow body trailing the deadweight depressor houses a CSS- 
developed low frequency synthetic aperture sonar and the 
SGMS. 

Reference 3 gives a detailed description of the feasi- 
bility demonstration test field layout, the testing procedures, 
and the performance of the acoustic and electro-optic sen- 
sors. Specific details relevant to the gradiometer demonstra- 
tion are given here. 


B. Gradiometer Results from 
the Feasibility Demonstration 


A linear target field was laid out to evaluate the SGMS 
performance. The field consisted of a row 200 m in length of 
eight small and medium-sized targets (mortar and artillery 
shells ranging in caliber from 60 up to 203 mm) running 
north-south; a second shorter row of three medium-sized tar- 
gets (two oil drums and a Mk82 500-lb bomb) parallel to, 
and 9 m east of, the first row; and a third row of two targets 
(a Mk83 1000-lb bomb and a Mk84 2000-lb bomb) 9 m east 
of the second row. A marker and a sonar calibration panel 
with ferrous anchors were laid at one end of the linear field. 
Estimates of the magnetic moments for the targets obtained 


from the D/C/L algorithm described in Section II.C are tabu- 
lated in Table I. These magnetic moments ranged from 0.03 
A-m? for one 60mm mortar shell up to 120 A-m? for the 
2000-Ib bomb. 

The detection range for these targets is also tabulated 
in Table I. Absolute range is given for a gradiometer with 
sensitivity of 3 pT/m-Hz"” at 0.1 Hz assuming a 10 dB sig- 
nal-to-noise ratio and relative range is given for fixed gradi- 
ent sensitivity (normalized to 1 for a 60mm shell with mag- 
netic moment of 0.03 A-m*). Observe that the 2000-Ilb bomb 
is detectable at 8 times the range of the 60mm mortar shell. 

The predicted results from the D/C/L algorithm to es- 
timate location and to classify the targets according to mag- 


TABLE | 
ESTIMATES OF MAGNETIC MOMENT FOR THE UXO TARGETS 





Target Moment Range 
(A-m*) (m) 
60mm shell 0.03 - 0.2 
81mm shell 0.3 
105mm sheil 0.7 
175mm shell 
203mm shell 
55 gal oil drum 
500-lb bomb 30 
1000-ib bomb 36 


2000- Ib bomb 48 


Relative 
Range 


1-1.5 





6-9 
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netic moment are displayed in Fig. 9 for one pass of the sys- 
tem through the target field. Open circles designate the ac- 
tual location of the targets, while the solid circles indicate 
the predicted positions. The size of the solid circles indi- 
cates the magnitude of the targets’ magnetic moment. For 
the feasibility demonstration, the nominal performance of 
the gradiometer channels was on the order of 3 pT/m-Hz'” 
at 0.1 Hz. The D/C/L algorithm was effective in localization 
and classification (by moment magnitude) all of the targets 
in this pass with the exception of the 1000-Ib Mk83 bomb. 
This exception provides an example of multi-target local- 
ization for which there are 2 targets (the 500-lb bomb and 
the 105mm shell) in this data window in addition to the 1000- 
lb bomb. In this case the algorithm successfully localized 
the 500-lb bomb and the 105mm shell with high signal 
strength, but failed to localize the more distant 1000-lb bomb 
with a relatively weak signal strength. The ferrous anchors 
for the marker and calibration panel at the right hand side of 
the map and two clutter objects in proximity to the 106mm 
shell (not a part of the target set) were also detected during 
this run. Double detections are displayed in Figure 9 for two 
targets, the 500- 1b bomb and one of the two oil drums. These 
represent target detections by the D/C/L algorithm in two 
separate data segments. The high degree of overlap for the 
double detections is suggestive of the degree of accuracy 
obtained with this algorithm. 


IV. Clutter Kejection 


Generally active acoustic approaches have proven to 
be an effective means to detect, classify and localize teth- 
ered sea mines or bottom mines proud with respect to the 
bottom in deeper waters. However the shallow-water bot- 
tom mine environment is an especially difficult acoustic en- 
vironment in which to operate. Interfering reverberations 
from the air/sea and sea/bottom interfaces, bottom topo- 
graphical features, general harbor debris, and mine burial 
present a difficult acoustic environment for bottom mine 
detection. In coastal regions, the density of debris clutter may 


TABLE I 
Features USED TO TRAIN THE NEURAL NeTwork DEVELOPED IN | 10! 





Magnetic Features Acoustic Features 
T 





Total magnetic moment Average signal-to-noise ratio 


(SNR) over object's length 

Average SNR over object's width 

Maximum of the ratio of SNR for 
length to SNR for width 


Y-component of moment 


Z-component of moment 


Depth esumaic 


Confidence level for correct object 
classification and localizanon 


| Number of intense pixels 
Target length 








| Target width 
| Estimated clutter density 
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Figure 10 

The probability of detection and classification for acoustic and 
magnetic sensors alone and for neural-network fusion data 
from both sensors given as a function of the number of false 
alarms per image in the two tracks from 10 to 27 meters on 
either side of the vehicle. 
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lead to a high false-alarm rate using conventional imaging 
sonar approaches alone. 

For effective clutter rejection, it is very desirable to 
use distinctly different sensor approaches. The application 
of two or more collocated sensors operating simultaneously 
has the potential to reduce false alarms and provide robust 
detection in a wide variety of background conditions. For 
mine reconnaissance and hunting, the combination of mag- 
netic sensors with sonars provides such an alternative. In the 
MADOM sea testing, more than 90% of the acoustically 
mine-like clutter was not magnetically mine-like. 

Several investigations have been conducted recently 
using automated neural network approaches to assess the 
merit of magnetic and acoustic data fusion [9], [10]. The 
following result was obtained courtesy of L. Smedley and 
G. Dobeck [10]. A set, consisting of 215 sonar images con- 
taining an assortment of drill targets and clutter objects, was 
assembled using data collected from sea tests with the SGMS 
and the MADOM low frequency synthetic aperture sonar. 
For each sonar image, the magnetic detections were co-reg- 
istered. An attractor-based k-nearest neighbor neural network 
was developed using the magnetic and acoustic features given 
in Table II. 

Receiver operating curves (ROCs) for this data were 
established for different detection ranges, recognizing the 
fact that the two sensors are effective over different ranges. 
The ROC for the individual sensors and fur the two-sensor 
data fusion (obtained from this neural network) are displayed 
in Fig. 10 for tracks on both sides of the vehicle in a range 
where both sensors were effective. For these tracks, there 
were 93 drill targets in the 215 images. 

We observe in this figure the improved detection and 
classification obtained using the data fusion. For an accept- 





able rate of 0.1 false alarms per image, the probability of 
detection and classification was 0.5 for the acoustic sensor 
alone and 0.77 for the magnetic sensor alone. The probabil- 
ity of detection and classification increased to 0.9 for a rate 
of 0.1 false alarms per image using this neural-network data 
fusion. Although these results were taken from a small data 
set of 215 images, we believe that the trend clearly demon- 
strates the substantial benefit of magnetic and acoustic data 
fusion for shallow water mine reconnaissance and UXO sur- 
vey. 


V. The Thin-Film 


Gradiometer 


We believe that the current technology, represented by 
the SGMS sensor, is reaching its performance limit. This 
technology is largely characterized by the use of bulk and 
wire niobium (Nb) superconducting components. Advances 
in Nb thin film technology to obtain increased low frequency 
sensitivity and the relative simplicity of the thin-film pro- 
cessing in contrast to labor-intensive assembly of bulk 
SQUID packages and the hand winding of wire loops are 
appealing. For mobile applications, the greater intrinsic bal- 
ance, i.e., common-mode rejection of the Earth’s magnetic 
field, obtained from thin-film lithography compared to 
manual winding of wire loops and the removal of bulk mag- 
netic components, including superconducting diamagnetic 
components such as shield canisters, is important to reduce 
anomalous signals in the gradiometer sense loops arising from 
acceleration-induced relative motion of parts. A project to 
develop a high sensitivity, all thin-film gradiometer sensor 
for mobile deployment is being pursued by the CSS, IBM 
Research [11], Ball Aerospace [12], Lockheed-Martin [13], 





Figure 11 
The TFG Cryogenic Probe Assembly. 








Quantum Magnetics [14], and the Naval Research Labora- 
tory [15] for demonstration in the Joint Countermine Ad- 
vanced Concept Technology Demonstration. 


A. Cryogenic and Koom- 
Temperature Electronics 


Major advances in Nb thin-film fabrication technol- 
ogy has led to the development, for the first time, of high 
quality low frequency SQUID based magnetic sensors uti- 
lizing Nb-AlO _-Nb tri-layer technology on a 5” scale. This 
work has led to the production of totally unshielded gradi- 
ometers which have been successfully demonstrated to op- 
erate in the Earth’s magnetic field. 

A cryogenic probe assembly for high performance in 
mobile operation has been manufactured with 3 tensor gra- 
diometer circuits mounted on a single-crystal silicon rod and 
mounted to the dewar neck plug (Fig. 11). The gradiometer 
circuits consist of 2 counterwound magnetometer loops, each 
3.8 cm square with a baseline of 5.3 cm, monolithically 
coupled to the SQUID. The precision lithography in con- 
junction with a configuration in which the sense loops, the 
SQUID washers, and their modulation and feedback coils 
are all pattered as gradiometers has provided extreme bal- 
ance in order to maintain full sensitivity in the presence of 
field changes on the order of 1000 nT. 

A benchtop version of high frequency flux-lock loop 
(FLL) feedback electronics with a modulation frequency of 
16 MHz, a factor of 15 to 30 times the current frequency 
available with commercial electronics, has been developed 
in order to assure specified signal-to-noise ratios required 
from the all thin-film gradiometer channels (using air-core 
thin-film output transformers in place of wire-wound fer- 
rous-core transformers) and to provide a high bandwidth for 
electromagnetic interference immunity [16]. 


B. The Advanced Liquid Helium 
Dewar 


A dewar prototype, referred to as the Advanced Liq- 
uid Helium Dewar, was developed to assure that the dewar 
would not limit sensor performance (Fig. 12). A flexible de- 
sign approach supported by detailed thermal, mechanical and 
field calculations was pursued. Stringent material selection 
and magnetic screening standards were established. The 
materials typically were chosen to be as magnetically clean 
as possible, with residual magnetizations 10 orders of mag- 
nitude smaller than that for soft steel. The magnetic gradient 
stability and the eddy current stability are 100 times better 
than the preceding state-of-the-art established with the 
SGMS. 

An exchange gas cooling approach was pursued in 
place of the convective cooling utilized in the SGMS dewar. 
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Figure 12 
The Advanced Liquid Helium Dewar. 








New innovations for thermal management have been imple- 
mented - anti-slosh baffles, thermal filters to isolate sensor 
area from bath temperature fluctuations, aluminized-mylar- 
blankets for radiation shielding custom etched to reduce the 
eddy currents, and thermal networks of 99.999% pure alu- 
minum wire for temperature uniformity. Temperature stabil- 
ity in motions typical of tow operations is on the order of 1 
uK/Hz" at 0.1 Hz a factor three orders of magnitude better 
than for the SGMS dewar [17]. 


C. Field Testing of the 
Integrated Sensor 


The integrated sensor prototype has been evaluated 
recently under field conditions. For this testing, the benchtop 
room temperature FLL electronics (which have a significant 
magnetic signature) are operated off a 15-m cable outside 
the test facility (remotely positioned some 13 m away from 
the sensor so as not to limit the performance of the sensor in 
motion). Stationary measurements have been conducted in 
this test setup in order to establish a baseline for the motion 
testing; i.e., to quantify any deterioration of performance in 
motion. In this configuration for the sensor, white noise on 
the order of 50 fT/m-Hz'” (10 u®/Hz'”) has been demon- 
strated. The knee for I/f noise occurred at approximately | 
Hz and the noise floor rose to approximately 200 fT/m-Hz'” 
(40 u® /Hz"*) at 0.1 Hz. The nominal balance of the gradi- 
ometers is estimated at | x10*/m (as measured indirectly from 
a comparison of compensated and uncompensated motion 
spectra obtained in preliminary motion testing), a factor of 
at least 20 greater than obtained with the wire-wound SGMS 
loops without the trimming procedures required for wire 
loops. The dynamic range of the gradiometers is 6x 10°. 

For comparison, these electronics evaluated in a labo- 
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ratory environment (integrated with a preliminary version 
of the gradiometer circuit currently being evaluated) operat- 
ing off a 2-m cable has white noise on the order of 4 w® /Hz'” 
with about 1/3 of the noise from the preamplifier. The knee 
for /f noise also occurred at approximately 1 Hz and the 
noise floor rose to approximately 7 u®,/Hz'” at 0.1 Hz. The 
electronics demonstrated a very high closed loop bandwidth 
exceeding 2.5 MHz and a very high slew rate greater than 
1x 10°@,/sec at frequencies up to 1 MHz [16]. In comparison, 
the SGMS has a bandwidth of 100 kHz and a slew rate of 
2.5x10°®,/sec. 


D. Development of a Field- 
Deployable TFG 


The existing TFG is being upgraded to a ruggedized 
field-deployable version. A fully-populated 5-channel cryo- 
genic probe assembly is being developed to replace the 3- 
channel laboratory unit. A compact field-deployable room- 
temperature electronics package is being developed to re- 
place the benchtop FLL electronics currently utilized for the 
laboratory version and to provide automated sensor control 
signal digitization, and data linking. Miniaturization of the 
electronics into a single integrated unit mounted onto the 
sensor (as required to obtain sensitivity of the integrated sen- 
sor in motion) represents a major undertaking. The entire 
analog and digital electronics for 5 gradiometer and 3 mag- 
netometer channels is being packaged into a unit 43 cm in 
diameter and 56 cm in length. The package must have mag- 
netic signature consistent with the sensitivity requirements 
in motion, power reduced by a factor of 30 compared to the 
laboratory prototype, and dimensions to minimize the length 
of the sensor body section for underwater deployment as 
much as possible. Production is in progress with final as- 
sembly and testing to be completed in 1997 for integration 
in an underwater tow system for mine hunting, specifically 
for demonstration in the Joint Countermine Advanced Con- 
cept Technology Demonstration, and for UXO surveys. 


VI. High T. 
Superconducting 
Technology 


As a result of nitrogen cooling, the development of 
sensors utilizing the high-T _ materials with nitrogen cooling 
provides an opportunity for significant size reduction, an ease 
of maintainability and convenience in comparison to the low- 
T, technology with helium cooling, factors critical to gain 
widespread acceptance of the superconducting technology 
over other magnetic sensor approaches. 





A. Perspective on Nitrogen 
Cooling for Naval Operations 


A broad-based assessment of refrigeration technology 
including liquid, solid, and triple-point nitrogen dewars and 
active cryocoolers was conducted [1], [18]. The conclusion 
from this assessment is that liquid nitrogen dewars represent 
the best choice for near-term development of a high perfor- 
mance high T. superconducting gradiometer for mobile ap- 
plications. The design for an open (vented) liquid nitrogen 
dewar with dimensions 45 to 75 cm in length and 30 to 50 
cm in diameter has been established consistent with sensi- 
tivity goals (Fig. 13). These dimensions are consistent with 
available space in underwater tow bodies of interest for U.S. 
Navy applications. A final choice for dimensions in any fi- 
nal dewar design would be based on a tradeoff between space 
and hold time. Results of the concept analysis indicate that a 
dewar with dimensions of 43 cm in diameter and 75 cm in 
length would have a hold time of approximately 33 days. 
This hold time is over 6 times greater than the hold time for 
the Advanced Liquid Helium Dewar which has the same di- 
ameter but is 150 cm long, twice the length of the nitrogen 
dewar. 

The benefits for naval mobile applications which can 
be obtained from these reduced cryogen requirements in- 
clude: (1) a significant reduction in down time during opera- 
tions; (2) affiliated reductions or elimination in labor require- 
ments for cryogen support during critical phases of an op- 
eration; (3) reduced failures in the cryogenic circuits or in 
the dewar (such failures typically occur during cryogen re- 
cycling); and (4) the elimination of an additional footprint 
on ship deck required for helium storage. In addition, the 





Figure 13 

Liquid nitrogen dewar concept with design versatile for 
multiple applications. This dewar has an outer diameter of 
43 cm (17”) (compatible with operation in a 53-cm (21”) 
vehicle) and a length of 75 cm (30”). The dewar is projected 
to have a hold time of 33 days. 
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use of liquid nitrogen in place of liquid helium significantly 
reduces supply logistics as a result of the wide availability 
of nitrogen on the market at domestic and most foreign ports 
and the availability of a large number of liquifiers in the U.S. 
Fleet, with at least 54 units identified. Significant cost sav- 
ings are expected from reduced costs for cryogen supply. 


B. Device Development under 
this Project 


Since 1993, there have been a number of laboratory 
results reported on magnetometer prototypes with white noise 
better than 200 fT/Hz'”. A number of test samples, magne- 
tometer circuits and gradiometer circuits have been devel- 
oped in conjunction with this project. This included a report 
of 26 fT/Hz'” at 1 Hz for a 2x2-cm magnetometer [19]. As 
an element of this program, the impact of flux trapping for 
unshielded operation of high T, sensors in the Earth’s mag- 
netic field has been investigated and identified to be more 
problematic than for the corresponding low T_ niobium thin- 
film sensors. Approaches are being pursued to circumvent 
current limitations in high T, fabrication technology [20]- 
[24]. In particular, the three-sensor gradiometer approach 
described in Section VII.B is an example of one means to 
circumvent these limitations. In that approach, magnetic-field 
coils are utilized to null out the Earth’s magnetic field at the 
sensing circuit. The field nulling significantly reduces noise 
associated with non-ideal magnetization effects in the high- 
T, superconducting material. 


C. 3-Axis High T. 
Magnetometer for Stationary 
Applications 


A 3-axis magnetometer prototype developed by 
Conductus under a SBIR contract has been evaluated (Fig. 
14). This sensor is being developed for stationary operation 
with a focus on geophysical applications. Results obtained 
to date are very promising. The sensor has been operated in 
the field totally unshielded without any deteriorated perfor- 
mance compared to its performance shielded. Magnetom- 
eter performance of 140 fT/Hz'” at 1 Hz has been demon- 
strated under stationary field conditions. This performance 
is a factor over 50 times better than the performance ob- 
tained from the best commercial fluxgate magnetometers. 
The 3 High-T, SQUID magnetometer circuits in this unit 
have been working reliably for over one year without failure 
[25]. 
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Figure 14 
3-channel high-T magnetometer. 








Vil. The Three-Sensor 
Gradiometer 


One advantage of the niobium-based superconducting 
technology, especially for mobile operation, is the ability to 
fabricate large scale counter-wound gradiometer sense loops 
using either niobium wire or multi-layer thin films with cross- 
overs. This allows signal subtraction using very low-noise 
passive circuits prior to signal processing with active ampli- 
fier circuits, which greatly reduces dynamic range require- 
ments from the active electronics. A gradiometer can be con- 
figured using two independent magnetometers with signal 





Figure 15 

The three-sensor gradiometer concept for one case in which 
there are two primary high T. SQUID-based magnetometers 
and a third fluxgate magnetometer for field nulling. 
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subtraction performed at the output of the two magnetom- 
eters. The CSS pursued this approach in the early 1970's 
using fluxgate technology. Good stationary performance was 
obtained at that time, but there was insufficient dynamic range 
in the processing of the differential signals to operate the 
sensor in motion. 

In order to circumvent this dynamic-range limitation 
arising by differencing two individual magnetometer signals 
during mobile operation, a novel approach patented by IBM 
Research is being pursued [26], [27]. A third magnetometer 
is used for common mode rejection, feeding back a signal to 
the two primary magnetometers which nulls out the ambient 
background field. This concept is denoted as the three-sen- 
sor gradiometer (TSG). The concept is depicted in Fig. 15 
for one case in which there are two primary high T. SQUID- 
based magnetometers and a third fluxgate magnetometer for 
field nulling. 


A. The Fluxgate Version of the 
TFG 


The basic TSG concept has been successfully pursued 
using room temperature fluxgate magnetometers (in place 
of the SQUID magnetometers). A laboratory prototype of 
the fluxgate TSG has been developed and demonstrated (Fig. 
16). This sensor features four commercial 3-axis triad sets 
of fluxgate magnetometers in a planar square array with a 
one-foot diagonal baseline. One triad set of fluxgate magne- 
tometers serves as a reference to measure the 3 mean mag- 
netic field vector components at the array. The remaining 
three sensor triads are each mounted inside their own 3-axis 
Helmholtz coil sets which null the mean magnetic field at 
the sensor triad. The residual signals at the triads are pro- 
cessed through their commercial electronics and appropri- 
ate combinations are subtracted via differential amplifiers. 
In this manner, 6 tensor gradient terms can be calculated of 
which 5 are independent. In the laboratory prototype, ana- 
log electronics are utilized to implement the magnetic feed- 
back currents with manual feedback adjustments determined 
to an accuracy of 5 decimals by a digital signal-processing 
routine. Sensitivity better than 0.3 nT/m at 0.1 Hz in motion 
has been demonstrated [28], [29]. 

A Phase I SBIR has recently been awarded to Quan- 
tum Magnetics to develop a ruggedized field-deployable 
version designed to improve the performance of the labora- 
tory prototype. This version will be compact and light weight 
for man-portable operation. It will feature an integrated com- 
puter for fully automated sensor control and signal process- 
ing and a display to provide the operator easy target detec- 
tion, classification, and localization. 

This sensor offers the opportunity to become the main- 
stay for man-portable magnetic surveys, replacing the total 
field magnetometer by offering unambiguous detection, 
moment classification, and localization. Although its range 








Figure 16 
Laboratory prototype of a 5-channel fluxgate gradiometer 
utilizing the three-sensor gradiometer concept. 








will be limited by factors of 4 to 10 compared to the super- 
conducting gradiometers described previously, this type of 
sensor will provide a low cost approach convenient for many 
applications and will avoid a need for cryogens. 


B. The High T_ Version of the 
TFG 


A project to develop and to demonstrate feasibility of 
a compact field-deployable high-T. superconducting gradi- 
ometer concept for mobile operation is being pursued by the 
CSS, IBM Research, Ball Aerospace, Quantum Magnetics, 





Figure 17 

Concept for the laboratory prototype of a 3-channel high-T. 
superconducting gradiometer utilizing the three-sensor 
gradiometer concept. 
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the Naval Research Laboratory, and Lockheed-Martin. As 
mentioned previously, it has been possible to fabricate high- 
sensitivity gradiometers in niobium technology by using 
monolithic wire or thin-film counterwound sense loops. 
Neither wire or thin-film monolithic loops can currently be 
manufactured using the high-T, technology to provide the 
higher sensitivities in motion offered by niobium technol- 
ogy. A high-T, gradiometer preliminary design has been es- 
tablished based on the TSG concept and its performance in 
motion has been modeled. The TSG approach circumvents 
the current limitations in high T, manufacturing technology, 
providing long baselines by using normal metal wire to con- 
nect the two SQUID magnetometers. 

Further improvement in performance is expected us- 
ing the high T. SQUID magnetometers in place of fluxgate 
magnetometers as a result of their intrinsically greater sensi- 
tivity. A laboratory test article is being developed to evalu- 
ate this concept (Fig. 17). This test article incorporates two 
flux cubes with each flux cube consisting of 3 orthogonal 
high-T, SQUID-based magnetometers with dimensions of | 
cm?. The two cubes, separated by a baseline of approximately 
30 cm, permits the synthesis of 3 tensor gradient compo- 
nents. For this test article, a 3-axis fluxgate triad is used to 
provide the 3 reference channels required in the feedback 
loop to null out the mean ambient background field. This 
high T, gradiometer is expected to surpass significantly the 
motion performance of any conventional non-superconduct- 
ing magnetic sensor technology and is projected to have sen- 
sitivity better than that of the low-T, SGMS. 


Conclusions and 
Summary 


Magnetic sensors provide one tool valuable for mo- 
bile search operations and surveys for targets with a signifi- 
cant magnetic signature. Superconducting SQUID-based 
sensors theoretically represent the most sensitive of known 
magnetic sensors. SQUID-based magnetometers have been 
demonstrated with sensitivities on the order of 1 fT/Hz'” at 
frequencies down to 0.1 Hz, while fluxgate and total-field 
magnetometers ale demonstrating sensitivities down to 1-10 
pT/Hz"” at 0.1 Hz. 

The U.S. Navy has developed the Superconducting 
Gradiometer/Magnetometer Sensor, a superconducting gra- 
diometer which has provided long-range detection compared 
to conventional non-superconducting magnetic sensors. This 
sensor has been utilized to demonstrate a capability for bur- 
ied mine detection and clutter rejection. As a result of the 
multi-channel approach, the sensor provides an accurate lo- 
calization capability and multi-target Giscrimination. The 
magnetic detection-and-classification signal processing de- 
veloped in conjunction with the sensor has proven to be ef- 
fective, providing a fully automated, real time capability. This 
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high-tech sensor has provided reliable, rugged performance 
in undersea tows conducted over a period of seven years. 
The sensor has been operated in the same tow vehicle adja- 
cent to a sonar without a loss in performance. This technol- 
ogy is available off-the-shelf to provide the greatest capabil- 
ity for magnetic detection and localization ever demonstrated. 
Work has continued recently with this sensor under the 
MUDSS Project to demonstrate its utility for UXO survey. 

A new approach incorporating all thin-film niobium 
components is being pursued for greater detection range in 
mobile operation. A laboratory prototype is being evaluated 
with white noise on the order of 50 fT/m-Hz"” rising to ap- 
proximately 200 fT/m-Hz"” at 0.1 Hz under stationary field 
conditions. A field deployable version is under development 
to be utilized in the Joint Countermine Advanced Concept 
Technology Demonstration in 1998. 

It is only a matter of time before high-T, devices will 
be available with performance comparable to that which has 
already been demonstrated with low-T_ devices. Magnetom- 
eter performance of 140 fT/Hz'” at 1 Hz has been demon- 
strated under stationary field conditions. This performance 
is a factor over 50 times better than the performance ob- 
tained from the best commercial fluxgate magnetometers. 
The benefits in reduced size, longer hold times, and reduced 
logistics and support make the high-T_ approach attractive 
compared to its low-T_ counterpart 

The localization capability afforded by 5-channel gra- 
diometers (and not previously afforded by conventional 
magnetic sensors) is expected to add impetus to the accep- 
tance of magnetics for mobile applications. The development 
of a 5-channel fluxgate gradiometer is currently in progress. 
Sensitivity better than 0.3 nT/m at 0.1 Hz in motion has been 
demonstrated. 

As a result of a much lower cost and no special sup- 
port requirements, such sensors will likely be sold in larger 
numbers for short-range applications. Their wider usage 


would then enhance the opportunity to display the utility of 


the greater classification and localization capability afforded 
by the 5-channel approach. With greater acceptance of the 
5-channel approach, end users will likely want to obtain sen- 
sors with greater sensitivity. It is likely that the eventual de- 
velopment of high-performance, reliable high-T. gradiom- 
eters will work synergisically with an increased acceptance 
of the more powerful magnetic signal-processing approaches. 
Hopefully, the time scales for this supply and demand will 
be commensurate and so expedite sensor developmental ef- 
forts. 
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Kesearch Notes 


The following are several of the major efforts the Navy is supporting which address the key priorities for mine 


countermeasures: 


Advanced Sensors for Unmanned Undersea and Airborne 
Vehicles 





Shallow Water 
Sensors 


Long Range 
Sensors 


Superconducting 


Synthetic Aperture Sonar Gradiometer 








Laser-iiuminated 
Electro-Optic 








Program develops technologies for remote reconnais- 
sance capabilities including high-resolution synthetic aper- 
ture sonar, electro-optical, and high critical temperature 
superconducting gradiometer sensors. 


Littoral Remote Sensing 

A program that will fuse the relevant data available 
from National sensors and standoff tactical sensor systems 
to provide the on-scene commander with geomorphologic/ 
environmental and mine/mine field location information. 


Shallow Water Airborne Mine Reconnaissance System 

A LIDAR system, a derivative of the "Magic Lantern" 
system, with improved bottom tracking, higher power and 
spatial resolution 


Alternative Detection, Classification, and Localization 
Development of new detection and classification ap- 
proaches using prototype neural network technology. 


Explosive Ordnance Disposal 

High resolution, hand held, imaging sonar systems for 
improved mine classification and identification capabilities, 
support to special warfare and explosive ordnance disposal. 


High Search Rate Sonar System 

A program to develop advanced sonar and signal pro- 
cessing technology for integration into unmanned undersea 
vehicles, including enhanced swath width toroidal volume 
search sonar, advanced side-looking sonar, and computer 
aided classification techniques. 


Advanced Degaussing Technology 

An effort aimed at significantly reducing steel-hull 
surface ship vulnerability to magnetic mines; focussing on 
closed loop systems and targeting new construction ships. 


Rapid Airborne Mine Clearance System Technology 
Demonstration 
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Develop and demonstrate the capability to rapidly tar- 
get and destroy mines in deep and shallow water from an 
airborne platform using supercavitating projectiles and tar- 
geting LIDAR. 
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Explosive Neutralization and Obstacle Clearance 


DET/SABRE Fire Control StS 





Thunder Road [Magic Carpet 


Efforts investigating technologies utilizing explosive 
nets, shaped charges, and explosively formed projectiles; de- 
ployment techniques, analytical performance prediction 
models and simulations are also developed. 


Improved Mechanical Mine Sweeping 


Advanced bottom contour following capability for 
sweeping close-moored mines 
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Environmental Support 

Critical to all aspects of mine countermeasures tech- 
nology development is an extensive knowledge of the envi- 
ronmental parameters that impact the performance of systems 
for detection, classification, and clearance of mines and ob- 
stacles. The efforts supported here include studies in atmo- 
sphere, ocean, ocean bottom and subbottom. Numerical 
predictive models and environmentally sensitive and adap- 
tive signal processing techniques are also supported. These 
lead to realistic simulation and training aids. 


Near-Term Mine Reconnaissance System 
A submarine launched unmanned undersea vehicle that 
will search for minelike objects and minefields clandestinely. 


Advanced Lightweight Influence Sweep System 

A system to provide broad and narrow band acoustic 
and magnetic energy sweep capability; using novel spark- 
gap and superconducting technology sources to emulate ship 
signatures. 











